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ABSTRACT 


A  parametric  sensitivity  analysis  of  "A  Regional  Input- 
Output  Model  For  Forecasting  Shore-Based  Navy  Workload"  in 
the  11th  Naval  District  is  presented  in  this  thesis.  The 
mathematics  of  the  Classical  Leontief  Input-Output  Model  is 
discussed.  The  formulation  of  "A  Regional  Input-Output  Model 
For  Forecasting  Shore-Based  Navy  Workload,"  by  the  Navy  Per¬ 
sonnel  Research  and  Development  Center,  San  Diego,  California, 
is  given.  The  model's  sensitivity  is  analyzed  by  changing  the 
input-output  coefficient  matrix  elements  in  each  row,  one 
row  at  a  time,  to  obtain  new  output  levels.  Resultant  per¬ 
centage  changes  in  the  output  levels  are  computed  and  recommen¬ 
dations  on  constructing  the  input  matrix  based  on  the  results 
of  this  analysis  are  discussed. 
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I.  INTRODUCTION 

In  1978,  the  Navy  Personnel  Research  and  Development 
Center,  located  at  San  Diego,  California,  formulated  "A 
Regional  Input-Output  Model  For  Forecasting  Shore-Based  Navy 
Workload"  in  the  11th  Naval  District.  At  the  time  of  its 
formulation,  this  model  was  designed  to  forecast  the  workload 
of  twenty  shore-based  activities  within  the  11th  Naval  Dis¬ 
trict.  This  thesis  is  the  result  of  a  parametric  sensitivity 
analysis  of  that  model.  In  the  sensitivity  analysis,  the 
author  examined  the  output  levels  of  the  model  when  the  input- 
output  coefficient  matrix  was  changed  by  a  small  amount  in 
each  row. 

Input-output  models  may  be  categorized  as  being  descrip¬ 
tive  or  prescriptive  in  nature  depending  upon  1)  the  goal  or 
purpose  of  the  model,  2)  the  applied  use  of  the  model,  3)  the 
structure  of  the  model,  and  4)  other  related  items.  Within 
this  broad  categorization,  it  is  useful  to  classify  models  in 
terms  of  their  use.  Input-output  models  may  be  classified  as 
models  of  pure  logical  consistency  or  models  of  pure  predic¬ 
tion.1  A  model  of  pure  logical  consistency  deals  only  with 
the  model's  logical  consistency.  It  addresses  the  question 
of  whether  or  not  certain  events  logically  lead  to  certain 


Jones,  C.R.,  A  Taxonomy  for  Naval  Force  Level  and  Structure 
Models ,  working  papers  Naval  Postgraduate  School,  Department 
of  Operations  Research  and  Administrative  Sciences,  Monterey* 
California,  April,  1971. 
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other  events.  There  are  no  empirical  tests  required  to  check 
this  class  of  model.  The  pure  prediction  model  generates 
forecasts  or  predictions  for  a  subsequent  time  period. 
Empirical  verification  of  this  class  of  models  may  be  done  by 
comparing  forecast  and  realized  observations.  "A  Regional 
Input-Output  Model  For  Forecasting  Shore-Based  Navy  Workload" 
in  the  11th  Naval  District  is  a  member  of  the  pure  prediction 
class  of  models. 

Section  II  of  this  thesis  provides  insight  into  how  an 
input-output  model  is  devised.  In  particular,  it  shows  how 
the  input-output  coefficients  are  computed  for  each  sector/ 
activity  of  a  general  input-output  model  and  gives  a  specific 
example  of  the  formulation  and  workings  of  a  very  simple  model 
The  mathematics  of  the  general  input-output  model  is  gone 
through  in  detail  and  the  model  is  reduced  in  form  to  a  simple 
expression  in  matrix  and  vector  notation. 

Presented  in  Section  III  is,  "A  Regional  Input-Output 
Model  For  Forecasting  Shore-Based  Navy  Workload, "  the  model 
formulated  by  the  Navy  Personnel  Research  and  Development 
Center.  A  step-by-step  breakdown  of  the  model's  formulation 
and  background  on  how  the  data  were  gathered  are  given.  Many 
years  of  data  collection  and  analysis  went  into  the  formula¬ 
tion  of  the  model.  This  section  presents  the  method  for  aggre¬ 
gation  of  sectors  and  coefficients. 

The  methodology,  chosen  by  this  author,  to  do  the  sensi¬ 
tivity  analysis  on  "A  Regional  Input-Output  Model  For  Fore¬ 
casting  Shore-Based  Navy  Workload"  in  the  11th  Naval  District 
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is  presented  in  Section  IV.  The  decision  to  do  a  parametric 
sensitivity  analysis  on  the  Input-Output  Coefficient  Matrix, 

A,  and  to  only  parametrize  rows  as  opposed  to  rows  and  columns 
or  some  combination  of  rows,  columns,  and  individual  elements 
of  the  input-output  coefficient  matrix  is  explained.  Also, 
the  choice  of  change  values,  for  parametrization,  of  between 
minus  25  percent  and  plus  25  percent  is  explained.  The  pro¬ 
cedure  for  computing  the  output  level  of  the  parametrized 
row/sector  is  outlined  as  is  the  procedure  for  computing  the 
Baseline  Output  Level  used  in  the  calculation  of  percentage 
changes  in  the  output  levels  of  the  twenty  sectors  of  the 
model.  The  crux  of  this  thesis  is  the  sensitivity  analysis 
of  the  output  level  of  a  sector  based  on  a  parametric  change 
to  that  sector's  row  of  the  input-output  coefficient  matrix 
and  the  manner  in  which  this  was  done  is  outlined  in  this 
section. 

The  results  of  the  methodology  outlined  in  Section  IV, 
when  applied  to  "A  Regional  Input-Output  Model  For  Forecasting 
Shore-Based  Navy  Workload,"  are  given  in  Section  V.  Upper  and 
lower  bounds  on  the  percentage  change  in  the  output  level  of 
a  parametrized  sector  are  computed  and  graphed  for  the  twelve 
sectors  this  author  felt  to  be  of  concern.  These  graphs  are 
presented  as  Figures  3  through  14  of  Section  V.  The  lower 
bounds  showed  a  linear  relationship  with  the  delta  values 
used  for  parametrization  to  compute  the  percentage  change  in 
the  output  level  for  a  particular  sector.  The  same  was  true 
for  the  computed  upper  bounds . 
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Conclusions  and  recommendations  are  given  in  Section  VI. 
The  fact  that  the  Public  Works  Center  was  the  most  sensitive 
sector  of  the  model  and  how  this  was  determined  are  presented. 
Inter-sectoral  relationships  are  discussed  and  the  recommen¬ 
dation  to  reevaluate  the  input-output  coefficient  matrix  in 
terms  of  the  zero  entries  is  made. 
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II.  THE  GENERAL  INPUT-OUTPUT  MODEL  AND  RELATED  MATHEMATICS 


l . 


The  classical  Leontief  Input-Output  Model  consists  of 
three  basic  elements.  A  square  matrix,  B,  is  used  to  repre¬ 
sent  the  industrial  classification  of  an  economy.  Each  row 
of  the  matrix,  B,  contains  the  amount  of  a  particular  indus¬ 
try's  final  product  that  is  used  as  input  for  producing  the 
final  product  of  industries  represented  in  the  columns  of 
the  matrix.  A  separate  row  vector  is  used  to  denote  the 
amount  of  each  primary  factor,  or  input  such  as  labor,  re¬ 
quired  by  each  industry  in  the  economy.  Finally,  a  separate 
column  vector  is  used  to  provide  the  amount  of  exogenous 
demand  that  exists  for  each  final  product  1 11; 13] . 

As  a  simplified  example,  consider  the  two- industry  economy 
in  which  industry  one  produces  water  and  industry  two  produces 
electricity.  If  b^  denotes  the  quantity  of  the  ith  industry’s 

i.U 

final  product  used  by  the  j  industry,  then  it  follows  that 
bn  represents  the  amount  of  water  used  in  the  production  of 
water  and  b12  represents  the  amount  of  water  used  in  the  pro¬ 
duction  of  electricity.  In  a  physical  sense,  b^  represents 
the  water  used  to  cool  and  supply  evaporators,  other  plant 
devices,  and  other  related  activities  of  the  two  producing 
industries.  In  a  likewise  manner,  b2^  and  b22  represent  the 
amount  of  electricity  used  to  operate  lights  and  other  elec¬ 
trical  components  in  the  water  and  electrical  industries 
respectively.  The  components  b^  b12,  b21'  and  b22  comPrise 
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the  square  matrix,  B.  Inasmuch  as  this  is  a  two- industry 
economy  example,  the  amount  of  labor  required  can  be  repre¬ 
sented  by  Xq^  and  and  the  exogenous  demand  for  each  final 
product  by  d^  and  d£.  The  total  real  output  of  the  two 
industries  can  be  represented  as  x^  and  X2 .  Figure  1  is 
that  of  the  general  model  for  an  n-industry  economy. 

The  details  of  solving  the  two-industry  economy  problem 
are  presented  as  an  introduction  to  the  methodology  used  in 
working  with  an  input-output  model .  The  solution  procedure 
deals  with  the  problem  of  determining  the  appropriate  levels 
of  output  at  which  the  two  industries  should  operate  such 
that  production  just  satisfies  the  total  demand  for  a  final 
product.  Simply  stated,  equilibrium  exists  in  the  two-industry 
economy  when  the  total  outputs  of  industry  one  and  those  of 
industry  two  are  in  balance  so  that  just  enough  of  each  is 
produced  to  satisfy  both  the  final  demand  and  the  input  re¬ 
quirements  for  each  product.  In  the  general  case,  the  two- 
industry  economy  becomes  an  n-industry  economy  and  the  solution 
is  analogous  to  the  two-industry  discussion. 

To  solve  the  two-industry  problem,  let  x^^  and  x2  represent 
the  total  outputs  of  industries  one  and  two  respectively. 

These  quantities  are  represented  as  row  sums  in  Figure  1.  It 
is  assumed  that  the  input-output  coefficients  are  constant 
[10]  and  that  the  transformation  of  inputs  into  outputs  is 
accomplished  in  a  linear  fashion.  From  the  discussion  of  the 
two- industry  economy  it  follows  that  for  each  unit  of  the  jth 

i.  u 

commodity  produced,  a  fixed  input  of  the  i  commodity  is 
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Figure  !•  The  Classic*!  Leontief  Model 
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required.  In  other  words,  a  unit  of  output  of  water  or 
electricity  requires  a  fixed  input  of  water  or  electricity 
respectively. 

The  input-output  matrix.  A,  may  be  derived  from  the  inter¬ 
industry  relationships.  In  general,  the  input-output  coeffi¬ 
cients  are  defined  by 


b .  ./x  . , 

ID  D 


where  a^  represents  the  amount  of  output  of  industry  i  used 
in  the  production  of  one  unit  of  output  of  industry  j,  b^^ 
represents  the  input  by  industry  j  to  industry  i,  and  x^ 
represents  the  total  real  output  by  industry  j.  Thus,  the 
input-output  coefficients  for  the  two- industry  problem  are 

all  *  bll/xl'  a12  "  b12/X2 '  a21  "  b21/Xl '  ***  a22  =  b22/x2’ 
The  total  output  requirements,  x^  and  X2  may  now  be 

expressed  as: 


a21Xl  +  a22X2  +  d2* 


Solving  these  two  equations  for  d^^  yields: 


1 


and 


S2  *  -a21Xl  +  »--a22)x2- 

Expressing  these  algebraic  expressions  in  matrix  notation  is 
as  follows: 
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The  above  matrix  notation  may  be  compactly  written  as 

(I-A)x  *  d. 

This  compactly  written  matrix  notation  will  be  used  through¬ 
out  the  thesis  with  A  denoting  the  matrix  of  a^s,  x  denoting 
the  vector  of  outputs,  5  denoting  the  vector  of  exogenous 
demands,  and  I  denoting  the  identity  matrix  (the  nxn  matrix 
with  Is  on  the  principal  diagonal  and  Os  elsewhere) .  Solving 
the  above  expression  for  x  completes  the  problem  and  it 
follows  that  x  =  (I-A)~*d  .  Now,  it  can  be  seen  that  the 
appropriate  level  of  output  for  each  of  the  two  industries 
is  uniquely  determined  for  a  given  exogenous  demand. 

The  assumptions  relating  to  the  Classical  Leontief  Input- 
Output  Model  will  now  be  addressed.  The  first  assumption 
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relating  to  a  single  output  for  each  sector  poses  no  concep¬ 
tual  problems.  However,  the  assumption  relating  to  the  con¬ 
stancy  of  input-output  coefficients  has  been  the  subject  of 
major  criticism  of  the  model.  A  very  enlightening  discussion 
relating  to  the  constancy  of  the  technological  coefficients  is 
provided  by  Hatanaka  [10] .  Hatanaka's  first  area  of  criticism 
relates  to  weaknesses  resulting  from  ignoring  certain  factors 
considered  in  production  theory.  He  specifically  categorizes 
shortcomings  as  being  the  result  of  ignoring  or  overlooking 
one  or  more  of  the  basic  factors  of  production.  First,  the 
model  ignores  the  occurrence  of  price  substitution.  Price 
substitution,  however,  is  a  very  real  phenomenon  since  choice 
among  alternative  inputs  is  very  dependent  on  relative  input 
prices.  Second,  no  accounting  is  made  for  economies  or  dis¬ 
economies  of  scale  for  inputs  which  are  consumed  by  output 
elements.  Third,  certain  factors  of  production,  which  are 
outside  the  model,  are  ignored.  This  criticism  results  from 
a  failure  of  the  model  to  account  for  such  things  as  capital 
stock  depletion  of  natural  resources.  Fourth,  the  model  fails 
to  account  for  joint  production.  This  shortcoming  results 
from  the  establishment  of  a  single  output  measure  for  each  of 
the  model's  input  elements.  An  example  of  this  joint  produc¬ 
tion  is  the  production  of  wool  and  mutton.  During  the  lifetime 
of  a  sheep,  the  hair  from  its  body  is  sheared  and  this  hair 
is  converted  into  wool.  As  the  sheep  matures  and  is  full 
grown,  it  is  slaughtered  and  the  flesh  of  the  animal  is  then 
offered  up  as  mutton.  Thus,  the  sheep  growing  industry  outputs 
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a  joint  product;  wool  from  young  sheep  and  mutton  from  more 
mature  sheep.  Fifth,  and  last,  the  model  fails  to  account 
for  technological  progress  during  the  time-frame  under  con¬ 
sideration  as  the  input-output  coefficients  are  assumed  to 
remain  constant. 

A  second  area  of  criticism  relates  to  weaknesses  that 
may  occur  as  a  result  of  the  manner  in  which  an  actual  model 
is  constructed.  Difficulties  may  arise  as  a  result  of  the 
manner  in  which  industries  are  defined  and  can  be  affected 
by  the  time  when  data  are  measured.  The  result  can  be  that 
the  input-output  coefficients  may  change  as  a  result  of  the 
manner  in  which  the  model  is  aggregated,  or  as  a  result  of  a 
change  in  the  mixture  of  old  and  new  production  procedures. 

To  visualize  the  effect  of  the  manner  of  aggregation,  assume 
that  there  are,  in  fact,  two  constant  coefficients.  If  these 
coefficients  are  added  during  aggregation  of  the  model,  there 
is  no  longer  any  reason  to  assume  that  the  constancy  of  indi¬ 
vidual  coefficients  has  been  preserved.  The  mixing  of  old  and 
new  production  procedures  can  be  clearly  illustrated  by  con¬ 
sidering  an  industrial  firm  that  is  faced  with  a  reduction 
in  demand  for  its  product.  Assume,  for  the  purpose  of  this 
example,  that  an  industry  does,  in  fact,  operate  an  old  and 
a  new  production  process  in  producing  its  output.  When  the 
industry's  sales  decrease,  it  will  undoubtedly  cut  back  on 
production  by  the  old  process.  As  business  improves,  and 
sales  begin  to  increase,  the  firm  will  probably  choose  to 
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expand  by  increasing  its  new  process  and  then  use  the  old 
process  as  capacity  utilization  increases.  The  end  result 
will  be  that  the  relative  proportions  of  old  and  new  produc¬ 
tion  lines  will  have  changed  and  hence#  the  input-output 
coefficients  will  have  changed.  A  final  problem  may  result 
from  the  methodology  used  in  building  the  model.  This  relates 
to  what  Hatanaka  [10]  calls  the  dilemma  of  aggregation.  This 
problem  relates  to  the  two-sided  nature  of  input-output  analy¬ 
sis,  or  the  fact  that  one  must  study  the  industry's  providing 
role  of  furnishing  its  output  as  input  to  an  intermediate 
or  final  output  industry  and  also  the  industry's  consuming 
role  of  using  supplied  inputs  in  the  production  of  outputs . 

The  problem  of  model  validation  can  present  additional 
difficulties.  There  are  essentially  two  approaches  that  may 
be  taken  to  determine  the  ability  of  a  model  to  provide  accurate 
forecasts.  The  difference  in  these  approaches  is  related  only 
to  the  time  period  for  which  the  predictions  are  made.  One 
method  of  prediction  is  to  construct  a  model  now  and  predict 
subsequent  outputs.  The  other  possibility  is  to  use  histori¬ 
cal  data,  available  in  the  economy,  to  build  a  model  and  pre¬ 
dict  outputs  for  which  a  realized  value  already  exists.  "A 
Regional  Input-Output  Model  For  Forecasting  Shore-Based  Navy 
Workload”  in  the  11th  Naval  District  was  validated  using  the 
latter  approach. 
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III.  THE  NAVY  PERSONNEL  RESEARCH  AND  DEVELOPMENT 
CENTER'S  REGIONAL  INPUT-OUTPUT  MODEL 


The  Navy  Personnel  Research  and  Development  Center, 
located  at  San  Diego,  California,  using  the  techniques  of 
input-output  modeling, formulated  "A  Regional  Input-Output 
Model  For  Forecasting  Shore-Based  Navy  Workload."  This  model 
was  a  quantitative  tool  to  forecast  workload  at  various  11th 
Naval  District  shore  activities  based  on  a  specific  cleet 
configuration.  The  input-output,  I/O,  model  was  developed  as 
a  means  by  which  I/O  analysis  could  capture  and  accurately 
measure  the  multitude  of  workload  demands  Navy  activities  place 
on  each  other.  It  could  also  help  in  determining  the  feasi¬ 
bility  and  cost  of  developing  a  Navy-wide  I/O  model  and  it 
could  be  used  to  explore  any  mathematical  and/or  computational 
problems  that  might  be  associated  with  the  formulated  model 
[7]  . 

In  developing  a  system  for  allocating  shore-based  ma’npower 
resources  in  the  11th  Naval  District,  major  emphasis  was  placed 
on  the  design  of  an  I/O  model  to  forecast  the  workload  of 
shore/support  activities  based  on  size  and  configuration  of 
the  fleet.  By  organizing  Navy  shore  activities  into  an  I/O 
matrix,  the  extent  to  which  each  activity  depended  on  every 
other  activity  to  produce  support  could  be  quantified  using 
historical  data.  I/O  analysis  was  seen  as  a  tool  to  determine, 
for  example,  the  impact  on  Long  Beach  Naval  Shipyard's  workload 
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of  introducing  an  additional  destroyer  into  the  overhaul 
schedule.  Also,  it  could  help  estimate  the  increased  work 
necessary  at  the  Naval  Supply  Center,  San  Diego,  California, 
to  service  more  shipyard  requisitions  in  support  of  the  over¬ 
haul  of  the  additional  destroyer.  Thus,  both  direct  and 
indirect  effects  of  the  fleet  were  captured  in  the  model. 

The  regional  I/O  model  was  formulated  as  a  standard  linear 
program  [7:9]  with  the  objective  being  to  minimize  total  activity 
output  while  satisfying  both  final  and  intermediate  demands 
for  that  output.  The  I/O  sectors  represented  11th  Naval  Dis¬ 
trict  shore-based  activities  with  final  demands  represented 

2 

by  the  ships  and  aircraft  of  the  fleet. 

The  regional  I/O  model  formulation  was: 

-  m 

Minimi zeJ  Z  -  \  x. 

i=l  1 


subject  to: 


(I  -  A)  X  >_  d 
X  >  0 


2 

Shore-based  activities  not  included  as  I/O  sectors  were 
treated  as  final  demands  on  the  system. 

3 

The  I-A  matrix  was  a  square  matrix  and  so  the  solution  to 
the  model  had  to  be  a  point  solution.  Thus,  the  objective 
function  served  no  purpose  except  to  formulate  the  model  as 
a  linear  programming  problem.  Hence,  the  objective  function 
could  be  dropped  since  the  "subject  to"  conditions  determine  the 
unique  solution  to  the  model. 
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where: 


m  =  number  of  I/O  sectors. 


x^  =  total  required  output  level  of  sector  i, 
d^  =  final  demand  on  sector  i. 


A  =  [a. .]  =  amount  of  input  required  by  sector  j  from 

sector  i  to  product  one  unit  of  output  from 
sector  j ,  and 


I  =  identity  matrix. 


Final  demand,  d^,  for  a  sector  was  treated  as  a  constant 
on  the  right-hand  side  of  the  linear  programming  formulation. 
Therefore,  calculation  of  final  demand  for  a  sector  by  fleet 
size,  mix,  location,  and  operating  tempo  was  done  outside  the 
LP  model  and  then  summed  up. 

The  formulation  of  the  final  demand  model  was: 


di  =  j=i  k=i  tii  nJkh  ^kh  riJkh 


where : 


n , 
jkh 


=  number  of  ships  of  type  j,  homeport  k, 
and  status  h,3 


fcjkh 


average  time  spent  in  status  h  by  a  ship 
of  type  j  and  homeport  k. 


4 

A  ship  would  be  in  one  of  the  following  statuses:  in 
port,  in  overhaul,  on  local  operations,  or  on  extended  cruise. 
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ri  jkh 


average  demand  rate  for  sector  i  goods  of 
ships  of  type  j,  homeport  k,  and  status  h. 


J 


number  of  ship  types. 


K 


number  of  homeports,  and 


H  =  number  of  statuses  (e.g.,  "in-port"  or 

"over-haul" ) . 


Nine  major  11th  Naval  District  shore  activities  were 
studied  in  detail  to  develop  an  I/O  model  representing  the 
fleet-support  demand  network.  The  activities  were  selected 
for  their  (1)  wide  range  of  functions,  outputs,  and  data 
problems,  (2)  manpower  intensities,  and  (3)  direct  and  indirect 
linkages  to  the  fleet.  An  I/O  formulation  of  the  fleet-support 
problem  required  the  creation  of  economic  sectors,  each  of 
which  produced  a  single,  unique  output.  While  the  workload 
of  many  shore  activities  could  be  characterized  by  a  single 
output  measure,  other  activities  produced  several  distinct 
outputs.  Thus,  it  was  necessary  to  divide  multioutput  activi¬ 
ties  into  several  sectors,  which  increased  the  number  of  I/O 
sectors  in  the  model  to  twenty.  Demands  placed  on  the  system 
by  units  not  included  in  this  set  of  twenty  sectors  are  all 
treated  as  final  demands  on  the  system.  This  includes  other 
11th  Naval  District  shore  activities.  Navy  activities  in 
other  geographical  areas,  and  the  operating  units,  ships  and 
aircraft,  of  the  fleet.  A  list  of  the  twenty  I/O  sectors  of 
the  model  and  the  units  that  were  used  to  measure  their  total 
workload  are  shown  in  Table  I. 
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Table  I 


I/O  Model  Sectors 


Sector/Activity  Output  Measure 

1  Naval  Supply  Center  Requisitions 

San  Diego 

2  Long  Beach  Naval  Shipyard  Man-days  of  Ship  Repair 


3  San  Diego-based  Intermediate 

Maintenance  Activities 

4  Naval  Air  Rework  Facility 

North  Island 

5  Public  Works  Center  (PWC) 

San  Diego 

6  PWC — Facilities  Maintenance 

7  PWC — Transportation 

8  PWC— Utilities 

9  Naval  Regional  Medical 

Center,  San  Diego 

10  Naval  Regional  Medical 

Center,  San  Diego 

11  Naval  Station — Military 

Personnel 

12  Naval  Station — Port  Services 

13  NAS  Miramar — Air  Operations 

14  NAS  Miramar — Aircraft 

Intermediate  Maintenance 

15  NAS  Miramar — Supply 

16  NTC — Recruit  Training 

17  NTC — Service  Schools 


Man-hours  of  Intermediate  • 
Maintenance 

Manhours  of  Air  Rework 

Man-hours  of  Direct  Labor 

Man-hours  of  Facilities 
Maintenance 

Gallons  of  Fuel 

Million  BTU's  of  Energy 

Admissions'* 

Outpatient  Visits'* 

Military  Personnel  Actions 

Port  Service  Hours 
Air  Operations 
Man-hours  of  AIM 

Requisitions 

Students 

Students 
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Table  I  (Cont’d) 


18 

Manpower — Active  Duty  Navy 

Personnel 

19 

Manpower — Dependents  of 

Personnel 

Active  Duty  Navy 

20 

Manpower — civilian 

Personnel 

°A  weighted  average  of  Admissions  and  Outpatient  visits 
was  used  for  the  output  measure. 
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To  develop  an  I/O  coefficient  matrix  for  the  model,  it 
was  necessary  to  study  each  activity  in  the  model.  This  series 
of  investigations  enabled  the  determination  of  each  sector's 
level  of  output  and  the  distribution  of  that  output  among  the 
other  sectors  of  the  model,  as  well  as  the  units  which  make 
up  final  demand  [1; 2; 3 ; 4 ; 5; 6; 8 ; 12; 14 ] .  At  each  activity, 
average  demand  rates  were  calculated  over  the  most  recent 
period  of  time  for  which  data  were  available,  generally 
FY 75-77 . 

In  addition  to  the  flow  of  goods  and  services  between  sec¬ 
tors  of  the  model,  the  analysis  of  activities  that  serve  the 
fleet  directly  required  the  estimation  of  final  demand  rates 
for  the  San  Diego-based  operating  units  of  the  Pacific  Fleet. 
This  involved  examining  the  level  of  demand  placed  on  an 
activity  by  each  of  the  approximately  120  ships  homeported  in 
San  Diego.  These  demand  rates  were  categorized  by  ship  type 
and  status.  Similar  analyses  were  performed  to  determine  the 
level  of  workload  required  to  service  each  type  of  aircraft  at 
the  Naval  Air  Rework  Facility.  Using  these  rates,  the  model 
was  able  to  compute  the  effect  on  final  demand  and  total  out¬ 
put  caused  by  changes  in  the  size  and  deployment  status  of 
operating  units  of  the  fleet. 

Since  workload  at  some  of  the  Navy  activities  could  only 
be  measured  in  terms  of  two  or  more  outputs,  special  estimation 
techniques  had  to  be  used  to  model  the  distribution  of  goods  and 
services  among  model  sectors  that  were  a  part  of  the  same  Navy 
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activity.  It  was  not  feasible  to  attempt  to  trace  each  input 
purchased  by  an  activity  to  the  production  department  from 
which  one  of  its  outputs  came;  hence,  certain  assumptions  were 
made  concerning  the  internal  use  of  inputs  within  these 
activities . 

The  matrix  shown  as  Figure  2  illustrates  the  economic 
interactions  among  sectors  of  the  model.  The  marked  cells 
indicated  the  location  of  nonzero  elements  in  the  I-A  matrix. 
Thus,  a  nonzero  entry  in  element  (i,j)  of  the  matrix  indicated 
that  the  output  of  sector  i  was  used  directly  in  the  production 
of  sector  j's  output.  As  an  example,  element  (1,2)  of  the 
matrix  indicated  an  interaction  between  the  Navy  Supply  Center 
and  Long  Beach  Navay  Shipyard.  This  meant  that  at  least  one 
requisition  was  processed  by  the  Navy  Supply  Center  for  Long 
Beach  Naval  Shipyard  and  that  a  change  in  the  shipyard  workload 
would  in  some  way  impact  the  workload  at  the  supply  center. 

In  an  effort  to  give  the  reader  a  working  knowledge  of 
the  input-output  coefficient  matrix,  the  first  row  of  A  and 
the  first  column  of  A  will  be  derived.  First,  the  input  matrix, 
B,  was  constructed  by  collecting  data  on  inputs  to  each  of  the 
twenty  sectors  of  the  model  from  each  sector  of  the  model.  The 
first  row  of  the  B  matrix  will  now  be  constructed.  Table  II 
is  a  table  of  the  first  row  of  the  input  matrix,  B.  During 
the  base  period  of  FY75-77,  The  Navy  Supply  Center,  NSC, 
generated  6131  requisitions  internally  which  had  to  be  processed 
by  NSC.  Thus,  the  input  to  sector  one,  NSC,  from  sector  one. 
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Figure  2-  lith  Naval  District  Region*!  I/O  Model  Impact  Matrix. 
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Table  II 


I 


Row  One  of  the  Input  Matrix,  B 


Sectors 

1-  5 

6131 

83387 

98472 

81634 

1090 

Sectors 

6-10 

0 

0 

0 

131589 

0 

Sectors 

11-15 

0 

6220 

64894 

0 

0 

Sectors 

16-20 

9069 

4282 

0 

0 

0 

Read  across  the  rows. 


NSC,  was  6131  requisitions  and  this  number  was  placed  in  the 
B  matrix  as  b^  Similarly,  Long  Beach  Naval  Shipyard,  sector 
two,  generated  83387  requisitions  which  NSC  processed.  This 
number  was  placed  in  the  B  matrix  as  b^^  2<  Entries  b^  4, 

and  b^  5  were  similarly  derived.  All  requisitions  for  the 
Public  Works  Center  were  generated  on  the  administrative  side 
of  the  house  of  the  Public  Works  Center,  hence.  Public  Works 
Center — Facilities  Maintenance,  Public  Works  Center — Trans¬ 
portation,  and  Public  Works  Center — Utilities  had  no  input  of 
requisitions.  Thus,  b1  g,  b.^  and  b.^  g  were  zero.  The  Naval 
Regional  Medical  Center  at  San  Diego  generated  131589  requisi¬ 
tions  which  were  processed  by  NSC.  Since  requisitions  were 
only  generated  in  the  administrative  side  of  the  house  of  Naval 
Regional  Medical  Center,  all  requisitions  were  put  into  the 
admissions  sector,  sector,  nine,  and  so  the  Naval  Regional 
Medical  Center — Outpatients  sector,  sector  ten,  had  no  input. 
All  requisitions  submitted  to  NSC  from  the  Naval  Station  at 
Saui  Diego,  6220  requisitions,  were  processed  as  Naval  Station 


Port  Services  requisitions  and  so  Naval  Station  Military 
Personnel,  sector  eleven,  had  no  input.  As  a  result  of  this, 
bi  became  zero  and  became  6220.  NAS  Miramar  sub¬ 

mitted  a  total  of  64894  requisitions  to  NSC  and  all  were  pro¬ 
cessed  as  NAS  Miramar — Air  Operations,  sector  thirteen,  requi¬ 
sitions.  Thus,  b^  ^2  became  64894.  Since  all  NAS  Miramar 
requisitions  were  processed  as  sector  thirteen  requisitions, 

NAS  Miramar — Aircraft  Intermediate  Maintenace,  sector  fourteen, 
and  NAS  Miramar — Supply,  sector  fifteen,  had  no  inputs  to  the 
matrix.  Consequently,  b^  and  b1  15  were  zero.  Naval  Train¬ 
ing  Center — Recruit  Training,  sector  sixteen,  submitted  9069 
requisitions  to  NSC.  Naval  Training  Center — Service  Schools 
Command,  sector  seventeen,  submitted  4282  requisitions  to  NSC. 
The  last  three  sectors  of  the  model  were  accounting  sectors  to 
track  personnel  only  and  so  no  requisitions  were  submitted  from 
these  sectors.  Thus,  Manpower — Active  Duty  Navy,  sector 
eighteen.  Manpower — Dependents  of  Active  Duty  Navy,  sector 
nineteen,  and  Manpower — Civilian,  sector  twenty,  had  no  inputs 
to  the  B  matrix  for  NSC,  sector  one.  This  completed  the  first 
row  of  the  B  matrix,  the  row  that  represented  the  Naval  Supply 
Center. 

The  first  column  of  the  B  matrix  was  derived  in  a  similar 
manner  to  the  first  row  of  the  matrix  and  Table  III  is  a  table 
of  the  first  column  of  the  input  matrix,  B. 

The  Naval  Supply  Center  input,  b^  ^  remained  unchanged. 

No  man-days  of  ship  repair  from  Long  Beach  Naval  Shipyard, 

sector  two,  were  used  by  NSC,  sector  one,  thus  b-  .  was  zero. 

z ,  ± 
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Table  III 


Column  One  of  the  Input  Matrix,  B 


Sectors  1-  5 

6131 

0 

0 

0 

0 

Sectors  6-10 

11734 

69983 

88906 

0 

0 

Sectors  11-15 

0 

0 

0 

0 

0 

Sectors  16-20 

Read  across  the 

0 

rows. 

0 

21 

0 

740 

Similarly,  no  man-hours  of  intermediate  maintenance  from  San 
Diego-based  Intermediate  Maintenance  Activities,  sector  three, 
were  used  by  NSC.  No  man-hours  of  air  rework  from  the  Naval 
Air  Rework  Facility  at  North  Island,  sector  four,  were  used 
by  NSC  and  no  man-hours  of  direct  labor  from  the  Public  Works 
Center,  sector  five,  were  used  by  NSC.  NSC  did  utilize  11734 
man-hours  of  facilities  maintenance  from  the  Public  Works 
Center — Facilities  Maintenance,  sector  six.  The  Public  Works 
Center — Transportation,  sector  seven,  contributed  69983  gallons 
of  fuel  to  the  operation  of  NSC.  In  order  to  heat  buildings 
of  NSC  and  run  auxiliary  equipments,  Public  Works  Center — 
Utilities,  sector  eight,  provided  88906  million  BTU's  of  energy 
to  NSC.  The  only  other  sectors  to  provide  inputs  to  NSC  were 
Manpower — Active  Duty  Navy,  Sector  eighteen,  21  personnel,  and 
Manpower — Civilian,  sector  twenty,  740  personnel.  Since  no 
other  sectors  had  inputs,  b^  ^ ,  b1Q  b^  1 ,  b^2  ^  b13 

^14  1'  ^15  1'  ^16  1'  ^17  1'  an<^  ^19  1  were  assigned  the 
value  zero.  This  completed  the  first  column  of  the  B  matrix. 
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The  remainder  of  the  B  matrix  was  filled  in  using  the  same 
techniques  as  just  shown  for  row  one  and  column  one. 

With  the  B  matrix  completed,  the  next  task  was  to  compute 
the  input-output  coefficients  and  put  them  in  the  A  matrix. 
This  was  done  by  dividing  the  inputs  in  the  columns  of  the  B 
matrix  by  the  total  real  output,  during  the  base  period,  of 
the  corresponding  sectors.  Table  IV  is  a  table  of  the  first 
row  of  the  input-output  coefficient  matrix,  A. 


Table  IV 

Row  One  of  the  Input-Output  Coefficient  Matrix,  A 


Sectors 

1-  5 

0.00580 

0.10500 

0.03410 

0.01260 

0.00560 

Sectors 

6-10 

0.0 

o 

• 

o 

0.0 

0.45510 

0.0 

Sectors 

11-15 

0.0 

0.62920 

0.59300 

0.0 

O 

o 

Sectors 

16-20 

1.64000 

0.90800 

o 

o 

0.0 

0.0 

Read  across  the 

rows . 

To  illustrate  the  procedure  described  above,  some  of  the 
elements  of  the  input-output  coefficients  for  the  Naval  Supply 
Center,  sector  one,  row  one  of  the  A  matrix  will  now  be  cal¬ 
culated.  To  compute  a.  ,  element  b..  .  was  divided  by  the  total  real 

1  #  1  i.  t  X 

output  for  NSC,  sector  one,  X^.  The  total  real  output  for 
each  of  the  twenty  sectors  may  be  found  in  Table  VI. 

ax  1  =  bx  l/Xl  =  6131/1057143  =  0.00580. 
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The  element  a.  _  was  computed  by  taking  the  element  b.  0  and 
dividing  it  by  the  total  real  output  for  Long  Beach  Naval 
Shipyard,  sector  two,  which  gave 

a^  2  =  2^2  =  ^3387/794165  =  0.10500  requisitions/man-day. 

The  element  a^  3  was  computed  by  taking  the  element  b.^  3  and 
dividing  it  by  the  total  real  output  for  San  Diego-based 
Intermediate  Maintenance  Activities,  sector  three.  This  gave 

a^^  3  =  b^  3/^3  =  98472/2887755  =  0.03410  requisitions/man-hour 

Since  the  elements  b16,  b1>7,  b1|8,  bl>10,  b1>n,  bl>14,  b1>15 

bl  18'  bl  19'  and  bi  20  were  2ero  elements  a^  a^  7, 

al , 8 '  al,10'  *1,11'  al, 14 '  al, 15 '  al , 18 '  *1, 19 '  and  al , 20 
were  assigned  the  value  zero.  This  completed  the  first  row 

of  the  input-output  coefficient  matrix,  A. 

To  calculate  the  values  of  the  input-output  coefficients 
in  the  first  column  of  A,  similar  logic  was  used.  Table  V  is 
a  table  of  the  first  column  of  the  input-output  coefficient 
matrix,  A.  The  first  element,  a^  ^  was  computed  above  and  so 
it  did  not  have  to  be  recomputed.  The  element  ag  ^  was  com¬ 
puted  by  taking  bg  ^  and  dividing  it  by  the  total  real  output 
for  NSC,  sector  one,  which  gave 

ag  ^  »  bg  ^/X^  *  11734/1057143  =  0.01110  man-hours/requisition 


Table  V 


Column  One  of  the  Input-Output  Coefficient  Matrix,  A 


Sectors 

1-  5 

0.00580 

O 

• 

0 

O 

O 

O 

• 

O 

0 

• 

0 

Sectors 

6-10 

0.0110 

0.6220 

0.08410 

0.0 

0.0 

Sectors 

11-15 

0.0 

0.0 

0.0 

0.0 

0.0 

Sectors 

16-20 

0.0 

0.0 

0.00002 

0.0 

0.00070 

Read  across  the  rows. 


The  element  a2g  ^  was  computed  by  taking  b2Q  ^  and  dividing 
it  by  the  total  real  output  for  NSC,  sector  one,  as 

a20  1  =  b20  ]/xi  =  740/1057143  =  0.00070  personnel/requisition 


The  elements  b2  ^ ,  b3  ^ ,  b4  ^ ,  b,.  ^ ,  bg  ^ ,  b1Q  ^ ,  b^1  ^ , 

b12 , 1 '  b13 , 1 '  b14 , 1 '  b15 , 1 '  b16 , 1'  b17 , 1 '  and  b19 , 1  were  a11 
zero.  Thus,  the  corresponding  a^j's  were  assigned  the  value 

zero . 

When  all  b. .  values  were  divided  by  X.  the  total  real 

13  3, 

output  of  sector  j,  the  complete  input-output  coefficient 
matrix.  A,  was  constructed.  When  the  input-output  coefficient 
matrix  was  subtracted  from  the  identity  matrix,  the  matrix  of 
technological  coefficients,  I-A,  was  generated.  The  techno¬ 
logical  coefficient  matrix  and  the  final  demand  vector  of  the 
model  are  provided  in  Appendix  [A] . 


The  study  which  generated  the  elements  of  the  B  matrix 
also  generated  the  total  real  output  vector  for  the  twenty 


Table  VI 


The  Total  Real  Output  Vector,  X 


Row 

Sector 

Output 

1 

NSC  -  SAN  DIEGO 

1057143.00 

REQUISITIONS 

2 

LBEACH  NAVAL  SHIPYARD 

794165.00 

MAN-DAYS 

3 

INTER.  MAINT.  ACTIVITY 

2887755.00 

MAN-HOURS 

4 

NAVAIR  REWORK  FACILITY 

6478878.00 

MAN-HOURS 

5 

PUBLIC  WORKS  CENTER 

194661.00 

MAN-HOURS 

6 

PWC  -  MAINTENANCE 

1074960.00 

MAN-HOURS 

7 

PWC  -  TRANSPORTATION 

609420.00 

GALLONS 

8 

PWC  -  UTILITIES  ENERGY 

3751626.00 

MILLION  BTUS 

9 

NAVREGMEDCEN  -  ADMISS . 

289142.00 

ADMISSIONS 

10 

NAVREGMEDCEN  -  OUTPAT. 

498157.00 

OUTPATIENTS 

11 

NAVSTA  -  MILITARY  PERS , 

45684.00 

PERSONNEL  ACTIONS 

12 

NAVSTA  -  PORT  SERVICES 

9886.00 

PORT  SERVICE  HOURS 

13 

NAS  MIRAMAR  -  AIR  OPS 

109434.00 

AIR  OPERATIONS 

14 

NAS  MIRAMAR  -  AIMD 

1499711.00 

MAN-HOURS 

15 

NAS  MIRAMAR  -  SUPPLY 

572336.00 

REQUISITIONS 

16 

NTC,  -  RECTRACEN 

5530.00 

STUDENTS 

17 

NTC  -  SERSCCLCCM 

4716.00 

STUDENTS 

18 

MANPOWER  -  ACTIVE  DUTY 

75904.00 

PERSONNEL 

19 

MANPOWER  -  DEPS  OF  ACDU 

17856.00 

PERSONNEL 

20 

MANPOWER  -  CIVILIAN 

28209.00 

PERSONNEL 
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sectors  of  the  model.  It  is  noteworthy  that  the  row  sums  of 
the  B  matrix  did  not  equal  the  corresponding  elements  of  the 
total  real  output  vector.  This  was  because  the  twenty  sectors 
of  the  model  were  those  of  selected  shore-based  activities  only 
and  the  total  real  output  vector  contained  requirements  for 
the  remaining  shore-based  activities  and  operating  fleet 
units  as  well. 

Using  the  technological  coefficient  matrix  described  above 
and  the  final  demand  vector,  the  solution  to  the  model  was 
found  by  solving 


(I  -  A)  X  =  d. 

This  yielded  a  baseline  output  level,  the  output  level  without 
parametric  changes  being  made  to  the  input-output  coefficient 
matrix.  This  was  a  twenty  element  vector  with  the  elements  of 
the  vector  being  the  corresponding  outputs  of  the  sectors  of 
the  model.  This  baseline  output  level  is  presented  as  Table 
VII. 

In  this  section  of  the  thesis,  the  reader  has  been  given 
the  formulation  of  "A  Regional  Input-Output  Model  For  Fore¬ 
casting  Shore-Based  Navy  Workload"  in  the  11  Naval  District 
and  its  component  parts.  The  elements  of  the  input  matrix, 

B,  and  their  origin  have  been  explained.  The  elements  of 
the  input-output  coefficient  matrix.  A,  and  their  computation 
have  been  explained.  The  technological  coefficient  matrix. 
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Table  VII 


Baseline  Output  Levels 


Row 

Sector 

Output  Level 

1 

NSC  -  SAN  DIEGO 

1572270.00 

Requisitions 

2 

LB EACH  NAVAL  SHIPYARD 

794165.00 

Man-Days 

3 

INTER.  MAINT.  ACTIVITIES 

2957636.00 

Man-Hours 

4 

NAVAIR  REWORK  FACILITY 

6478878.00 

Man-Hours 

5 

PUBLIC  WORKS  CENTER 

2589904.00 

Man-Hours 

6 

PWC — MAINTENANCE 

1573271.00 

Man-Hours 

7 

PWC — TRANSPORTATION 

1503492.00 

Gallons  of  Fuel 

8 

PWC — UTILITIES  ENERGY 

5939721.00 

Million  BUT's 

9 

NAVREGMEDCEN — ADMISSIONS 

310200.00 

Admissions 

10 

NAVREGMEDCEN— OUTPATIENTS 

1726666.00 

Outpatients 

11 

NAVSTA— MILITARY  PERSONNEL 
ACTIONS 

45684.00 

Personnel  Actions 

12 

NAVSTA— -PORT  SERVICES 

9886.00 

Port  Service  Hours 

13 

NAS  MIRAMAR — AIR  OPERATIONS 

109434.00 

Air  Operations 

14 

NAS  MIRAMAR— AIMD 

1602596.00 

Man-Hours 

15 

NAS  MIRAMAR — SUPPLY 

620330.00 

Requisitions 

16 

NTC — RECRUIT  TRAINING 

CENTER 

5530.00 

Students 

17 

NTC — SERVICE  SCHOOLS 

COMMAND 

4716.00 

Students 

18 

MANPOWER— ACTIVE  DUTY  NAVY 

106825.00 

Personnel 

19 

MANPOWER — DEPENDENTS  OF 
ACTIVE  DUTY  NAVY 

140704.00 

Personnel 

20 

MANPOWER— C IV IL  IANS 

56891.00 

Personnel 

I-A,  has  been  presented  and  the  solution  to  the  model,  the 
baseline  output  level,  has  been  presented. 
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IV.  METHODOLOGY  USED  IN  THE  SENSITIVITY  ANALYSIS 


The  Navy  Personnel  Research  and  Development  Center's  "A 
Regional  Input-Output  Model  For  Forecasting  Shore-Based  Navy 
Workload"  in  the  11th  Naval  District  was  composed  of  a  twenty 
by  twenty  technological  coefficient  matrix,  I-A,  and  a  twenty 
element  final  demand  vector,  d.  In  order  to  perform  a  sensi¬ 
tivity  analysis  on  this  model,  it  was  decided  to  parametrically 
change  each  of  the  sixty-three  nonzero  a^’s,  i  =  1,  2,  . .., 

20;  j  =1,  2,  ...,  20,  of  the  input-output  coefficient  matrix, 
A,  by  a  small  percentage  called  delta  and  to  compute  the  output 
of  the  model  with  the  new  technological  coefficient  matrix, 
(I-A*),6  and  the  given  final  demand  vector,  d. 

The  elements  of  the  final  demand  vector  were  held  constant 
during  this  process  because  the  interest  was  in  what  would 
happen  to  the  output  of  the  model  if  the  input-output  coeffi¬ 
cients  were  parametrically  changed  for  each  row  of  A  and  a 
solution  was  obtained  to  the  model  with  the  final  demand  vector 
being  held  constant. 

The  decision  was  made  to  change  an  entire  row  rather  than 
to  change  individual  elements  of  the  input-output  coefficient 
matrix.  This  was  because  of  the  manner  in  which  the  a^’s  were 
computed,  as  discussed  in  Sections  II  and  III  of  this  thesis. 
The  input-output  coefficients  were  computed  from  the  input 


6The  asterisk  indicates  a  parametric  change  has  been  effected. 


matrix,  B.  The  elements  of  B  were  taker  from  studies  of  each 
sector  of  the  model,  and  thus,  an  entire  row  of  B  was  the 
result  of  one  study.  Any  one  column  of  B  was  the  partial 
result  of  twenty  such  studies.  It  was  more  likely  that  if 
over-counting  or  under-counting  was  done  in  the  sorting  of 
a  sector's  data  to  get  the  input  for  each  sector,  the  differ¬ 
ence  would  be  more  prevalent  in  the  row  than  in  the  column 
corresponding  to  the  sector  under  study.  Since  the  input 
matrix,  B,  reflected  this  difference  in  its  rows,  it  then 
followed  that  the  a^'s,  computed  from  the  b^'s,  would  preserve 
this  property.  From  the  relationship  between  a^^  and  b^, 
given  in  Section  I,  it  could  be  seen  that  parametrizing  a  row 
of  the  input-output  coefficient  matrix.  A,  by  a  percent  was 
equivalent  to  parametrizing  the  same  row  of  the  input  matrix, 

B,  by  a  percent.  That  is,  since 

a^j  -  b^j/total  real  output  of  industry  j, 
if 


a 


9 


then 


aij 


=  a  a.. 


as  long  as  the  total  real  output  of  industry  j  is  not  changed. 
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Further,  it  was  found,  on  inspection  and  testing,  that  para- 
metrization  of  rows  produced  consistently  larger  changes  than 
did  parametrization  of  columns.  Thus,  it  was  decided  to  con¬ 
centrate  on  the  effect  that  parametric  changes  to  the  rows 
would  make  and  not  to  concentrate  on  the  effect  that  parametric 
changes  to  the  columns  would  make. 

In  Section  III,  Table  VII,  the  baseline  output  level,  Q, 
of  the  model  was  presented.  Q  was  a  twenty  element  vector 
with  each  element  of  the  vector  representing  one  of  the  twenty 
sectors  of  the  Navy  Personnel  Research  and  Development  Center's 
model.  To  compute  output  levels  of  the  model  with  parametric 
changes  effected  to  the  A  matrix,  each  row  of  the  input-output 
coefficient  matrix  was  parametrically  changed  by  a  percentage, 
one  row  at  a  time,  and  a  new  output  level,  Q1,  was  obtained. 

One  Q1  for  each  row  of  the  input-output  coefficient  matrix 
was  calculated  for  each  delta.  The  elements  of  Q1  represented 
output  levels  for  each  sector  of  the  parametrized  model. 

To  calculate  the  percentage  change  in  output  level  for 
each  of  the  twenty  sectors,  the  baseline  output  level  of  that 
sector  was  subtracted  from  that  sector’s  element  in  Q1,  the 
difference  divided  by  the  baseline  output  level  of  that  sector, 
and  the  result  multiplied  by  100.  The  following  algebraic 
expression  was  utilized  to  compute  the  percentage  changes: 


percentage  change 
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lU 

in  the  output  of  the  j  sector  of  the  model  when  the  i 
row  of  the  input-output  coefficient  matrix  had  been  para¬ 
metrized;  i  =  1,  2,  ...  ,  20;  j  =  1,  2,  ...  ,  20.  Using  this 
expression,  twenty  percentage  changes  were  computed  for  each 
parametrized  row. 

It  was  necessary  to  decide  which  of  these,  if  any,  were 
significant  enough  to  use  for  later  analysis.  To  do  this,  a 
test-statistic  had  to  be  decided  upon.  It  was  noted  that  most 
percentage  changes  tended  to  be  at  or  near  zero.  It  was  de¬ 
cided  to  use  the  mean  percentage  change  as  the  test-statistic. 
The  mean  percentage  change  for  the  output  of  the  model,  given 
a  parametrized  row,  was  computed.  This  mean  percentage  change 
was  considered  significant  for  the  delta  value  used  to  compute 
the  output  of  theparametrized  model.  Each  percentage  change, 
of  the  twenty  computed  above,  was  then  compared  to  this  test- 
statistic  and  if  it  was  greater  than  the  test-statistic,  it 
was  considered  significant.  The  significant  percentage  changes 
were  identified  and  used  to  further  inspect  the  inter-sectorial 
relationships  in  the  next  section  of  this  thesis.  If  the  per¬ 
centage  change  was  less  than  the  test-statistic,  it  was  con¬ 
sidered  not  significant  and  disregarded. 

As  an  example  of  this  methodology,  significant  percentage 
changes  in  the  output  level  of  the  model  with  row  nine,  the 
row  corresponding  to  Naval  Regional  Medical  Center — Admissions, 
being  changed  by  ten  percent  will  now  be  given.  Table  VIII 
is  the  vector  of  percentage  changes  for  the  model  when  row  nine 
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Table  VIII 


Vector  of  Percentage  Changes  in  Outputs  for  a 
Parametric  Change  in  Row  Nine  of  Ten  Percent 


Sectors 

1-  5 

0.06 

o 

• 

o 

O 

• 

o 

O 

o 

0.03 

Sectors 

6-10 

0.02 

0.02 

0.06 

0.69 

0.11 

Sectors 

11-15 

0.0 

0.0 

0.0 

0.0 

0.0 

Sectors 

16-20 

0.0 

0.0 

0.16 

0.14 

0.13 

Read  across  the  rows;  values  in  percent. 

The  mean  percentage  change  equals  0.07  percent. 


of  the  input-output  coefficient  matrix  was  parametrized  by 

9 

ten  percent.  This  vector  was  written  as  Q . ,  j  =  l,  2,  ...  , 

20. 

The  significant  percentage  changes  are  those  greater  than 
the  mean  percentage  change.  Sequentially  going  through  the 
table,  the  significant  percentage  changes  are:  sector  nine, 
0.69  percent;  sector  ten,  0.11  percent;  sector  eighteen,  0.16 
percent;  sector  nineteen,  0.14  percent;  and  sector  twenty, 

0.13  percent. 

Interpretation  of  these  percentage  changes  is  now  easily 
done.  If  the  actual  input-output  coefficients  for  the  Naval 
Regional  Medical  Center — Admissions,  sector  nine,  were  ten 
percent  larger  than  the  estimates,  one  could  expect  the  output 
levels  for  each  of  the  sectors  to  be  the  given  percentages 
greater  than  the  baseline  output  levels.  For  example,  the 
output  for  Manpower — Civilian,  sector  twenty,  would  be  56964 


i 

i 
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or  0.13  percent  greater  than  the  56891  personnel  indicated 
by  the  baseline  output  level. 

The  methodology  employed  in  this  sensitivity  analysis 
was  such  that  a  minus  delta  implied  an  over-estimation  of 
the  input-output  coefficients  by  the  given  percentages  and  a 
positive  delta  implied  an  under-estimation  of  the  input-output 
coefficients.  It  was  reasonable  to  assume  that,  in  the  worst 

7 

case,  delta  should  be  in  the  range  -10  percent  to  +10  percent. 
Thus,  the  computed  percentage  changes  given  in  Section  V  of 
this  thesis  were  for  various  values  of  delta  in  the  range  -25 
percent  to  +25  percent.  This  wide  range  encompassed  the  reason¬ 
able  range  and  allowed  for  a  wider  range  over  which  to  do  the 
analysis . 

Depending  by  how  much  a  particular  row  of  the  input-output 
coefficient  matrix  was  likely  to  be  over-estimated  or  under¬ 
estimated,  lower  or  upper  bounds  on  the  percentage  change  in 
the  output  level  of  the  Navy  Personnel  Research  and  Development 
Center's  model  were  established.  Decreasing  the  coefficients 
of  a  row  by  ten  percent  and  solving  the  model  would  yield  an 
output  level  for  each  of  the  twenty  sectors.  Similarly,  in¬ 
creasing  the  coefficients  of  that  same  row  by  ten  percent  would 
yield  another  output  level  for  each  sector.  The  percentage 
changes  computed  for  each  sector  would  then  be  bounds  on  the 


7 

These  values  were  arrived  at  in  consultation  with  Mr. 
Thomas  Blanco  of  the  Navy  Personnel  Research  and  Development 
Center  who  did  much  of  the  estimation  of  the  input-output 
coefficients. 
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percentage  changes  of  the  respective  sectors.  As  an  example, 
the  coefficients  of  row  nine,  the  row  corresponding  to  the 
Naval  Regional  Medical  Center — Admissions,  were  decreased  by 
ten  percent  and  percentage  changes  in  the  output  levels  of  the 
twenty  sectors  were  calculated.  The  percentage  change  for 
sector  nine  was  -0.69  percent.  When  the  coefficients  were 
increased  by  ten  percent  the  percentage  change  for  sector  nine 
was  0.69  percent.  These  two  percentage  changes  were  lower 
and  upper  bounds  on  the  percentage  change  for  sector  nine  for 
any  values  used  for  parametrization  between  plus  and  minus 
ten  percent.  Decreasing  the  input-output  coefficients  of  row 
nine  by  five  percent  caused  a  percentage  change  of  -0.34 
percent  and  increasing  the  coefficients  of  row  nine  caused  a 
percentage  change  of  0.34  percent.  Both  of  these  percentage 
changes  were  bounded  by  -0.69  percent  and  0.69  percent.  There¬ 
fore,  as  long  as  the  change  in  the  input-output  coefficients 
for  row  nine  was  within  plus  and  minus  ten  percent,  the  resul¬ 
tant  percentage  change  in  the  output  of  sector  nine  would  be 
bound  by  plus  and  minus  0.69  percent.  Further,  the  bounds 
computed  for  the  parametrized  sector  were  also  bounds  on  the 
percentage  change  for  all  other  sectors  of  the  model. 
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V.  RESULTS  OF  THE  PARAMETRIC  SENSITIVITY  ANALYSIS 


In  order  to  conduct  a  parametric  sensitivity  analysis,  a 
computer  algorithm  was  written  following  the  methodology  out¬ 
lined  above  in  Section  IV.  First,  the  algorithm  computed  the 
baseline  output  level,  the  solution  to  the  Navy  Personnel  Re¬ 
search  and  Development  Center's  Regional  Input-Output  Model 
without  parametrically  changing  the  rows  of  the  input-output 
coefficient  matrix.  Next,  each  row  of  the  input-output  coeffi¬ 
cient  matrix  was  parametrically  changed  by  a  percentage  and  a 
new  output  level  was  calculated.  This  gave  twenty  new  output 
levels,  one  for  each  of  the  twenty  sectors  of  the  model.  Delta 
values  of  -25  percent  to  +25  percent  were  used  in  +5  percent 
increments.  Thus,  a  total  of  ten  values  were  used  for  para¬ 
metrically  changing  each  row  of  the  input-output  coefficient 
matrix.  This  provided  ten  output  levels  for  each  sector  to 
be  compared  with  the  baseline  output  level  of  that  sector . 

The  comparison  was  done  by  computing  the  percentage  change 
from  the  sectors  baseline  output  level  and  identifying  the 
significant  percentage  changes. 

Certain  sectors  of  the  model,  by  the  definition  of  signifi¬ 
cant  percentage  change  given  in  Section  IV,  caused  more  signi¬ 
ficant  percentage  changes  than  others.  Since  the  requirement 
for  manpower  was  a  primary  concern  of  the  Navy  Personnel 
Research  and  Development  Center's  model,  particular  attention 
was  given  to  sectors  which  when  their  input-output  coefficient 


row  was  changed  caused  significant  percentage  changes  and  to 
all  manpower  sectors.  This  group  included  Naval  Supply  Center — 
San  Diego,  sector  one;  San  Diego-based  Intermediate  Maintenance 
Activities,  sector  three;  Public  Works  Center,  sector  five; 

PWC — Maintenance,  sector  six;  PWC--Transportation,  sector 
seven;  PWC — Utilities,  sector  eight;  Naval  Regional  Medical 
Center--Admissions,  sector  nine;  Naval  Regional  Medical  Center-- 
Outpatients,  sector  ten;  NAS  Miramar — AIMD ,  sector  fourteen; 

NAS  Miramar — Supply,  sector  fifteen;  Manpower — Active  Duty 
Navy,  sector  eighteen;  Manpower--Dependents  of  Active  Duty 
Navy,  sector  nineteen;  and  Manpower--Civilian,  sector  twenty. 

A  summary  of  the  significant  percentage  changes  caused  by 
parametrizing  a  row  of  the  input-output  coefficient  matrix 
corresponding  to  a  particular  sector  is  given  as  Table  IX. 

The  sectors  that  had  no  significant  percentage  changes  were 
ones  whose  row  of  the  input-output  coefficient  matrix  consisted 
of  all  zeros.  This  meant  that  these  sectors  had  no  input 
from  the  other  sectors  of  the  model.  Thus,  their  rows  of 
the  input  matrix  were  all  zero.  Since  the  input  matrix  rows 
and  the  input-output  coefficient  matrix  rows  were  all  zero, 
on  parametrization  the  coefficients  remained  zero.  Solution 
of  the  parametrized  model  then  gave  the  same  output  levels  as 
the  baseline  output  level  and  thus  no  significant  percentage 
changes  were  computed.  The  outputs  of  these  sectors  were 
always  constant. 

Line  numbering  on  the  figures  of  this  section  was  done  to 
indicate  the  sector  the  line  corresponded  to.  The  lines  were 
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Table  IX 


Summary  of  Significant  Percentage  Changes  Caused  by 
Parametrizing  a  Particular  Sector 


Sector  Parametrized 

Sectors  with  Significant 
Percentage  Changes 

1 

NAVAL  SUPPLY  CENTER 

1,  7 

2 

LBEACH  NAVAL  SHIPYARD 

None 

3 

INTER.  MAINT.  ACTIVITIES 

1,  3 

4 

NAVAIR  REWORK  FACILITY 

None 

5 

PUBLIC  WORKS  CENTER 

5,  7 

6 

PWC— MAINTENANCE 

5,  6,  7 

7 

PWC — TRANSPORTATION 

5,  7 

8 

PWC — UTILITIES  ENERGY 

5,  8 

9 

N AVREGME  DC  EN — ADMISSIONS 

9,  10,  18,  19,  20 

10 

NAVREGMEDCEN — OU  TP AT I ENTS 

10 

11 

NAVSTA— MILITARY  PERSONNEL 

None 

12 

NAVSTA — PORT  SERVICES 

None 

13 

NAS  MIRAMAR — AIR  OPS 

None 

14 

NAS  MIRAMAR — AIMD 

14,  15 

15 

NAS  MIRAMAR— SUPPLY 

15 

16 

NTC — RECTRACEN 

None 

17 

NTC — SERSCOLCOM 

None 

18 

MANPOWER— ACTIVE  DUTY  NAVY 

10,  18,  19 

19 

MANPOWER — DEPENDENTS  OF  ACDU 

10,  19 

20 

MANPOWER — C IV IL IAN 

20 

provided  so  the  reader  could  have  a  visual  representation  of 
the  percentage  change  computed  for  a  known  value  by  which  the 
input-output  coefficient  matrix  row  for  the  sector  in  the  figure 
heading  was  changed.  Further,  to  visualize  the  bounds,  the 
reader  may  select  a  percentage  change  for  a  sector  of  concern. 


The  bound  will  be  the  ordinate  of  the  least-squares  line  at 
the  point  on  the  line  corresponding  to  the  chosen  percentage 
change . 

The  least-squares  lines  on  some  figures  were  very  flat, 
indicating  small  slopes.  This  meant  that  the  percentage 
changes  resulting  from  changes  in  the  values  of  the  rows  of 
the  input-output  coefficient  matrix  were  very  small.  Some 
of  these  "flat"  lines  may  seem  uninteresting  to  the  reader 
but  they  were  included  for  the  many  managers  who  may  read 
this  thesis.  A  quick  look  at  a  graph  is  sometimes  easier  to 
understand  than  to  have  to  compute  values  and  do  the  graphing 
for  one's  self. 

2 

The  coefficient  of  determination,  r  ,  was  computed  to 
establish  how  well  the  delta  values  and  their  associated 

percentage  changes  for  a  significant  sector  fit  the  linear 

2  2 
regression.  AT  r  equal  zero,  you  have  no  fit  and  at  r  equal 

one,  you  have  a  perfect  fit.  In  the  case  of  all  sectors  which 

2 

showed  a  significant  percentage  change,  r  was  approximately 

2 

equal  to  one.  The  r  values  were  in  all  cases  greater  than 
0.9995.  Thus,  the  percentage  changes  computed  for  the  given 
delta  values  were  very  close  to  being  linear.  Since  the  model 
was  linear  and  the  parametric  changes  formed  linear  transfor¬ 
mations  of  the  input-output  coefficient  rows,  the  resultant 
percentage  changes,  computed  from  the  sectorial  outputs,  would 
themselves  be  expected  to  be  linear.  It  was  felt  that  the 
reason  the  regression  line  did  not  pass  exactly  through  the 
origin  was  because  of  round-off  by  the  computer  in  calculating 
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the  inverse  of  the  twenty  technological  matrices,  I-A* .  Some 
of  the  percentage  changes  were  not  considered  to  be  significant 
and  so  the  regression  lines  for  their  sectors  were  not  graphed. 
Nonetheless,  the  reader  might  be  interested  in  the  magnitude 
of  the  nonsignificant  percentage  changes  so  all  percentage 
changes  resulting  from  parametrization  of  a  model  sector  are 
given  in  Appendix  [B] .  The  reader  may  check  this  appendix 
for  any  sector  of  concern. 

A.  NAVAL  SUPPLY  CENTER 

In  terms  of  significant  percentage  changes,  for  all  delta 
values  used  in  the  sensitivity  analysis,  parametric  changes 
to  the  input-output  coefficients  of  sector  one.  Naval  Supply 
Center — San  Diego,  produced  significant  changes  in  the  output 
levels  of  Navi  Supply  Center,  sector  one,  and  PWC — Transporta¬ 
tion,  sector  seven.  This  could  be  reasoned  out  since  in  order 
for  the  Naval  Supply  Center  to  serve  the  other  sectors,  supply 
them  with  goods,  the  Naval  Supply  Center  must  have  transporta¬ 
tion  to  transport  goods  from  the  supply  depot  to  the  requisi¬ 
tioning  activity.  Thus,  the  significant  change  in  the  output 
level  of  PWC--Transportation ,  given  the  change  in  the  Naval 
Supply  Center's  output  was  expected.  Other  sectors  showed 
change  as  well,  but  based  on  the  definition  of  "significant 
percentage  change,"  the  sectors  of  concern  were  the  Naval 
Supply  Center  and  PWC — Transportation.  The  significant  percen¬ 
tage  changes  that  occurred  as  a  result  of  parametrizing  row 
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one  of  the  input-output  coefficient  matrix  are  graphed  as 
Figure  3. 

B.  INTERMEDIATE  MAINTENANCE  ACTIVITIES 

Parametric  changes  to  the  nonzero  elements  of  the  input- 
output  coefficients  in  row  three,  the  row  corresponding  to 
San  Diego-based  Intermediate  Maintenance  Activities,  produced 
significant  changes  in  the  output  levels  of  Naval  Supply 
Center — San  Diego,  sector  one,  and  San  Diego-based  Intermediate 
Maintenance  Activities.  Knowing  that  the  Intermediate  Main¬ 
tenance  Activities  were  heavily  dependent  on  the  Naval  Supply 
Center  to  supply  repair  parts  for  maintenance,  this  linkage 
was  to  be  expected.  The  significant  percentage  changes  that 
resulted  from  parametrizing  row  three  are  graphed  as  Figure  4. 

C.  PUBLIC  WORKS  CENTER 

Changing  the  nonzero  elements  of  row  five  of  the  input- 
output  coefficient  matrix,  the  row  representing  the  Public 
Works  Center,  sector  five,  produced  significant  changes  in 
the  output  levels  of  Public  Works  Center,  sector  five,  and 
PWC — Transportation,  sector  seven.  Most  of  the  output  for  the 
Public  Works  Center  required  transportation  of  some  sort  to 
expend  the  man-hours  so  the  change  in  the  output  level  of 
PWC — Transportation  was  to  be  expected.  One  would  also  expect 
changes  in  the  output  levels  of  the  sectors  which  supported 
the  Public  Works  Center,  i.e..  Naval  Supply  Center--San  Diego, 
PWC — Utilities,  Manpower — Civilian,  and  as  expected,  changes 
occurred.  However,  by  the  definition  of  "significant,"  these 
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Least-Squares  Lines 

1  -  NSC — San  Diego 
Slopes  0*002  Intercepts  0-000 

3  -  Intermediate  Maintenance  Activities 
Slopes  0-024  Intercepts  0-000 


Delta  Values  In  Percent 


Figure  4-  Intermediate  Maintenance  Activities,  Sector  3- 


54 


changes  were  not  significant.  It  was  also  noted  that  the 
percentage  change  in  output  of  the  Public  Works  Center  was 
approximately  proportional  to  the  amount  of  change  made  to 
the  input-output  coefficients  of  that  sector.  Figure  5  is 
a  graph  of  the  significant  percentage  changes  resulting  from 
parametrizing  row  five. 

D .  PWC — MAINTENANCE 

Parametrically  changing  the  nonzero  elements  of  row  six 
of  the  input-output  coefficient  matrix,  the  row  corresponding 
to  PWC — Maintenance,  sector  six,  produced  significant  changes 
in  the  output  levels  of  Public  Works  Center,  sector  five, 

PWC- -Main ter ance ,  sector  six,  and  PWC — Transportation,  sector 
seven.  Since  the  model  was  originally  formulated  with  all  PWC 
sectors  combined  into  the  one  sector.  Public  Works  Center, 
the  change  in  the  output  level  of  the  Public  Works  Center  was 
expected.  Since  maintenance  was  required  to  keep  the  fiscal 
plants  of  the  Public  Works  Center  up,  the  change  in  the  output 
level  of  sector  six,  PWC — Maintenance,  was  expected.  The 
change  in  PWC — Transportation  was  expected  also  since  trans¬ 
portation  was  required  to  transport  personnel  and  equipment 
to  maintenance  tasks  in  the  11th  Naval  District.  Figure  6  is 
the  graph  of  significant  percentage  changes  caused  by  para¬ 
metrizing  row  six. 

E.  PWC — TRANSPORTATION 

Changing  the  elements  of  row  seven,  the  row  corresponding 
to  sector  seven,  PWC — Transportation,  evidenced  significant 
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Least-Squares  Lines 

5  -  Public  Works  Center 

Slope:  0-203  Intercept:  0-004 

6  -  PWC — Maintenance 

Slope:  0-320  Intercept:  0-00? 

7  -  PWC— Transportation 

Slope:  0-046  Intercept:  0-000 
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changes  in  the  output  levels  of  Public  Works  Center,  sector 
five,  and  PWC — Transportation,  sector  seven.  Inasmuch  as 
the  relationship  between  all  sectors  of  the  Public  Works 
Center,  as  explained  above,  existed,  these  changes  were 
expected.  The  graph  of  the  significant  percentage  changes 
resulting  from  parametrizing  row  seven  is  presented  as  Figure  7. 

F .  PWC — UTILITIES 

Changing  row  eight,  the  row  corresponding  to  PWC — Utilities, 
sector  eight,  caused  significant  changes  to  Public  Works 
Center,  sector  five,  and  PWC — Utilities,  sector  eight.  As 
explained  above,  the  relationship  of  all  Public  Works  Center 
sectors  would  lead  one  to  expect  these  changes.  The  graph  of 
these  significant  percentage  changes  is  given  as  Figure  8. 

G.  NAVREGMEDCEN— ADMISSIONS 

The  sector  whose  input-output  coefficient  row,  when 
changed,  caused  the  greatest  number  of  significant  percentage 
changes  was  the  Naval  Regional  Medical  Center — Admissions, 
sector  nine.  Changing  the  nonzero  elements  of  row  nine  pro¬ 
duced  significant  changes  in  the  output  levels  of  Naval  Regional 
Medical  Center--Admissions ,  Naval  Regional  Medical  Center — 
Outpatients,  Manpower--Active  Duty  Navy,  Manpower — Dependents 
of  Active  Duty  Navy,  and  Manpower--Civilian .  The  change  in 
the  output  level  of  Naval  Regional  Medical  Center--Admissions 
was  to  be  expected  since  the  input-output  coefficients  in  that 
row  were  changed.  A  change  in  the  number  of  personnel  admitted 
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to  the  Naval  Regional  Medical  Center  would  generally  cause 
a  change  in  the  number  of  outpatient  visits  and  so  a  change 
in  outpatient  visits  was  to  be  expected.  Given  that  active 
duty  personnel  and  dependents  of  active  duty  personnel  com¬ 
prise  a  large  part  of  the  number  of  patients  who  are  seen  at 
the  hospital,  change  in  the  output  of  the  sectors  corresponding 
to  these  two  groups  of  people,  sectors  eighteen  and  nineteen 
respectively,  was  expected.  Civilians  make  up  a  large  part 
of  the  staff  at  the  Naval  Regional  Medical  Center  and  so 
change  in  the  workload  at  the  hospital  would  precipitate  a 
change  in  the  civilian  workforce  at  the  hospital.  These  signi¬ 
ficant  percentage  changes  are  given  as  the  graph  of  Figure  9. 

H .  NAVREGMEDCEN— OUTPATIENTS 

Changing  row  ten  of  the  input-output  coefficient  matrix, 
the  row  which  corresponded  to  Naval  Regional  Medical  Center — 
Outpatients,  sector  ten,  showed  only  a  change  in  the  output 
of  that  sector.  The  Naval  Regional  Medical  Center  was  originally 
only  one  sector  in  the  model  and  so  it  was  assumed  that  change 
in  the  output  level  of  the  related  sector  was  absorbed  in  the 
Naval  Regional  Medical  Center — Admissions  sector,  sector  nine. 
Figure  10  is  the  graph  of  significant  percentage  changes. 

I.  NAS  MIRAMAR- -AIMD 

Changing  row  fourteen,  the  row  corresponding  to  NAS 
Miramar — AIMD,  sector  fourteen,  produced  significant  changes 
in  the  outputs  of  NAS  Miramar — AIMD,  sector  fourteen,  and  NAS 
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1 

Miramar- -Supply ,  sector  fifteen.  This  was  to  be  expected 
since  the  Aircraft  Intermediate  Maintenance  Depot  is  heavily 
dependent  on  NAS  Miramar--Supply  for  repair  parts  to  effect 
repairs  to  aircraft.  See  Figure  11  for  a  graph  of  the  signi¬ 
ficant  percentage  changes. 

J.  NAS  MIRAMAR — SUPPLY 

Changing  row  fifteen,  the  row  corresponding  to  NAS  Miramar — 

Supply,  sector  fifteen,  produced  significant  change  in  the 
output  level  of  NAS  Miramar — Supply,  sector  fifteen,  only. 

This  was  to  be  expected  since  all  requisitions  from  NAS 
Miramar  were  processed  by  NSC  as  NAS  Miramar — Air  Operations, 
requisitions.  The  percentage  changes  ranged  from  -1.94  percent 
for  a  delta  value  of  -25  percent  to  +1.94  percent  for  a  delta 
value  of  +25  percent.  Since  the  sector  under  study  was  the 
only  sector  to  show  significant  percentage  change,  no  graph 
is  provided. 

K.  MANPOWER — ACTIVE  DUTY  NAVY 

Parametrically  changing  the  nonzero  elements  of  row  eighteen, 
the  row  corresponding  to  Manpower--Active  Duty  Navy,  sector 
eighteen,  produced  significant  percentage  changes  in  the  output 
levels  of  Naval  Regional  Medical  Center--Outpatients ,  sector 
ten.  Manpower — Active  Duty  Navy,  sector  eighteen,  and  Man¬ 
power — Dependents  of  Active  Duty  Navy,  sector  nineteen.  Changing 
the  number  of  active  duty  personnel  in  the  11th  Naval  Dis¬ 
trict  should  cause  a  change  in  the  number  of  dependents  of 
active  duty  personnel.  With  a  change  in  the  number  of  active 
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duty  personnel  and  a  change  in  the  number  of  dependents  of 
active  duty  personnel  in  the  11th  Naval  District  one  would 
expect  a  change  to  occur  in  the  number  of  outpatient  visits 
at  the  Naval  REgional  Medical  Center.  Significant  percentage 
changes  are  graphed  as  Figure  12. 

L.  MANPOWER— DEPENDENTS  OF  ACTIVE  DUTY 

Changing  row  nineteen,  the  row  corresponding  to  Manpower — 
Dependents  of  Active  Duty  Navy,  sector  nineteen,  produced 
significant  changes  in  the  output  levels  of  Naval  Regional 
Medical  Center — Outpatients,  sector  ten,  and  Manpower-Dependents 
of  Active  Duty  Navy,  sector  nineteen.  This  was  to  be  expected 
since  dependents  of  active  duty  personnel  comprise  a  large 
part  of  the  Naval  Regional  Medical  Center’s  outpatient  service. 
See  Figure  13  for  the  graph  of  significant  percentage  changes. 

M.  MANPOWER— CIVILIAN 

Changing  row  twenty  of  the  input-output  coefficient  matrix, 
the  row  corresponding  to  Manpower — Civilian,  sector  twenty, 
produced  a  significant  change  in  the  output  level  of  Manpower 
Civilian,  sector  twenty.  One  would  expect  changes  to  occur 
in  the  output  levels  of  all  sectors  where  civilian  manpower 
was  utilized  but  the  model  only  showed  a  change  of  any  kind 
in  the  output  for  Manpower — Civilian.  Figure  14  is  the  graph 
of  this  significant  percentage  change. 

As  discussed  above,  upper  and  lower  bounds  may  be  taken 
from  the  graphs.  Choose  the  plus  and  minus  percentages  a  row 
is  most  likely  to  be  off  by.  Then  construct  the  necessary 
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perpendiculars  from  the  chosen  values  to  the  least-squares 
line  of  the  sector  of  interest.  From  the  points  of  inter¬ 
section  on  the  least-squares  line  to  the  ordinate  axis  of 
the  graph  construct  other  perpendiculars.  Where  the  normal 
lines  intersect  the  ordinate  axis,  bounds  may  be  read. 
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VI .  SUMMARY 


A.  CONCLUSIONS 

This  section  presents  the  conclusions  resulting  from  the 
parametric  sensitivity  analysis. 

The  most  sensitive  sector  to  parametric  change  of  its 
input-output  coefficients  was  the  Public  Works  Center.  When 
its  input-output  coefficient  row,  row  five,  was  parametrized 
and  percentage  changes  were  computed  for  the  various  values 
of  delta,  this  sector  had  a  least-squares  line  with  a  slope 
of  0.98.  This  implied  that  if  row  five,  the  row  of  the  A 
matrix  corresponding  to  the  Public  Works  Center  was  parametrized, 
the  percentage  change  in  the  output  of  that  sector  would  be 
approximately  equal  to  the  percentage  used  for  parametriza- 
tion.  Coupled  with  this  was  the  percentage  change  of  PWC — 
Transportation.  The  slope  of  the  least-squares  line  for  PWC — 
Transportation,  sector  seven,  was  0.23.  This  implied  that  if 
the  row  of  the  input-output  coefficient  matrix  representing 
the  Public  Works  center  were  changed  by  a  small  amount  then 
the  output  level  of  PWC — Transportation  would  be  approximately 
that  of  the  baseline  output  level  of  PWC — Transportation  plus 
this  baseline  output  level  times  0.23  times  delta.  This  inter¬ 
sectorial  relationship  indicated  that  the  managers  of  these 
two  sectors  would  want  to  be  very  mindful  of  changes  in  the 
output  level  of  the  Public  Works  Center  as  a  change  in  its 
output  level  would  precipitate  a  significant  percentage 
change  in  the  output  level  of  PWC — Transportation. 
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The  second  most  sensitive  sector  was  Manpower — Dependents 
of  Active  Duty,  sector  nineteen.  This  sector  had  a  least- 
squares  line  with  a  slope  of  0.88.  Thus,  changing  its  input- 
output  coefficient  row  would  cause  the  output  level  of  the 
model,  with  the  parametrized  row  to  be  approximately  that  of 
the  baseline  output  level  for  sector  nineteen  plus  the  baseline 
output  level  times  0.88  times  delta.  The  intersectorial 
relationship  between  this  sector  and  Naval  Regional  Medical 
Center — Outpatients,  sector  ten,  was  such  that  changing  row 
nineteen  of  the  A  matrix  produced  significant  percentage 
changes  in  the  output  level  of  Naval  Regional  Medical  Center — 
Outpatients.  Further,  changing  the  coefficients  representing 
Manpower- -Dependents  of  Active  Duty  produced  a  Naval  Regional 
Medical  Center — Outpatients  least-squares  line  with  a  slope 
of  0.27.  Hence,  changing  row  nineteen  would  cause  the  output 
of  sector  ten  to  be  approximately  that  of  the  baseline  output 
level  for  sector  ten  plus  the  baseline  output  level  times 
0.27  times  delta.  Thus,  a  change  in  the  number  of  dependents 
would  precipitate  a  change  in  the  number  of  outpatients  who 
visited  the  Naval  Regional  Medical  Center. 

The  third  most  sensitive  sector  was  Naval  Regional  Medi¬ 
cal  Center — Outpatients,  sector  ten.  For  this  sector,  a 
least-squares  line  slope  of  0.71  was  computed.  This  implied 
that  if  the  input-output  coefficients  of  row  ten  were  changed, 
the  output  of  the  model  with  the  changed  coefficients  would 
be  approximately  equal  to  the  baseline  output  level  for  sector 
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ten  plus  the  baseline  output  level  for  sector  ten  times  0.71 
times  delta. 

The  fourth  most  sensitive  sector  was  PWC — Transportation, 
sector  seven.  Changing  the  input-output  coefficients  of  row 
seven  yielded  a  line  with  a  slope  of  0.62  when  the  output 
level  of  the  model,  with  the  changed  coefficients,  was  com¬ 
pared  to  the  baseline  output  level  of  sector  seven.  Thus, 
changing  the  coefficients  of  row  seven  would  mean  that  the 
output  of  sector  seven  would  be  approximately  equal  to  the 
baseline  output  level  of  sector  seven  plus  the  baseline  out¬ 
put  level  of  sector  seven  times  0.62  times  delta.  Coupled 
with  this  was  a  significant  percentage  change  in  the  output 
level  of  Public  Works  Center,  sector  five.  Making  the  afore¬ 
mentioned  changes  to  row  seven  would  produce  a  least-squares 
line  for  the  Public  Works  Center  with  a  slope  of  0.19.  There¬ 
fore,  by  changing  row  seven  the  output  of  the  Public  Works 
Center  would  be  approximately  equal  to  the  baseline  output 
level  of  the  Public  Works  Center  plus  that  baseline  output 
level  times  0.19  times  delta.  In  terms  of  number  of  man-hours, 
this  number  was  not  too  great,  however,  by  the  definition  of 
significant  percentage  change,  it  was  significant. 

The  fifth  most  sensitive  sector  to  parametrization  was 
Manpower — Civilian,  sector  twenty.  Changing  the  coefficients 
in  row  twenty  produced  a  line  with  a  slope  of  0.50  when  the 
output  level  of  the  model,  with  the  parametrized  coefficients, 
for  sector  twenty  was  compared  to  the  baseline  output  level 
for  sector  twenty.  Thus,  changing  row  twenty  by  an  amount 
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delta  would  produce  an  output  level  approximately  equal  to 
the  baseline  output  level  for  Manpower — Civilian,  sector 
twenty,  plus  this  baseline  output  level  times  0.50  times 
delta. 

There  were  other  sectors  one  could  concentrate  on,  how¬ 
ever,  all  others  had  least-squares  line  slopes  of  less  than 
0.50.  The  reader  might  check  the  others  by  going  through  the 
figures  of  Section  V  and  comparing  the  slopes  of  the  least- 
squares  lines.  The  measure  of  sensitivity  used  to  determine 
the  most  sensitive  sectors  was  the  slope  of  the  least-squares 
line  through  the  points  generated  by  computing  percentage 
changes  at  the  given  values  of  delta. 

The  relationship  of  a  change  in  a  row  of  the  input-output 

coefficient  matrix  to  a  change  in  the  original  input  matrix 

was  a  simple  one.  The  input-output  coefficients,  a..'s,  were 

^-J 

computed  from  the  input  matrix,  B,  and  the  total  real  output 
vector,  X.  The  algebraic  expression  for  this  transformation 
was 


aij 


=  b . ./x . 
JO  1 


and  since  elements  of  the  total  real  output  vector  were  held 
constant,  it  was  simple  to  compute  the  inverse  of  this  trans¬ 
formation  given  a  change  to  a^.  To  show  the  inverse  transfor¬ 
mation,  assume  row  i  of  the  input-output  coefficient  matrix 
was  changed  by  an  amount  k.  Then 
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(l+k)ai>  =  b±  __  */x  j 

and  solving  this  expression  for  b^  yielded 


b.  *  =  (1  +  k)  a .  x.  . 

i  ■  i  •  j 


Now,  let's  assume  the  elements  of  row  one  of  the  A  matrix 

★ 

were  decreased  by  ten  percent  and  compute  new  values  of  b^  ^ 

★ 

where  b,  ,  would  be  the  value  of  b,  ,  had  the  original  b, 

1/1  1/1  1* 

been  such  that  the  values  of  a^  *  were  computed  from  them. 
Then, 


(l-0.10)a1#1  =  b1/1/x1 

which  on  substitution  would  be 


( 1  -  0 . 10)  ( 0 . 00580)  =  b1  ./1057143  requisitions 

■L  » 


or 


(0.90) (0.00580)  =  b.  ./1057143  requisitions. 

■L  ,  ■*■ 


Solving  this  would  yield 


b.  .  =  5518  requisitions 

1/1 
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Since  5518  is  90  percent  of  6131,  the  original  b 
shows  that  a  change  of  k  to  the  a^ _  implied  a  change  of  k 
to  the  corresponding  b^..  Simply  stated,  changing  a^  by 
some  amount  yielded  a  corresponding  change  of  the  same  amount 
to  the  respective  b.^  . 

The  results  of  this  analysis  may  be  used  to  evaluate  a 
sector's  output  based  on  a  historical  change  in  the  sector's 
performance.  Suppose  a  sector's  actual  inputs  from  the  other 
sectors  were  increased  by  some  percentage.  The  model's  output, 
when  the  model  was  changed  to  reflect  the  new  inputs,  could  be 
predicted  by  finding  the  amount  of  the  change  on  the  appropriate 
figure  or  in  Appendix  [B]  and  reading  the  percentage  change 
from  the  baseline  output  level . 

Suppose  managers  of  the  Naval  Supply  Center  and  the  Naval 
Regional  Medical  Center  were  to  ask  themselves  "what  if  the 
model's  input-output  coefficients  were  off  by  plus  and/or  minus 
ten  percent?  What  range  of  requisitions  or  admissions  might 
the  respective  sectors  encounter."  Thus,  the  sectors  managers 
would  be  interested  in  total  variability  of  forecasts  with 
the  given  model  input-output  coefficients. 

If  the  inputs  to  the  Naval  Supply  Center  were  actually  ten 
percent  greater  in  magnitude  than  the  inputs  in  the  original 
model,  the  percentage  change  in  the  output  of  the  Naval  Supply 
Center  would  be  3.30  percent.  This  would  imply  that  instead 
of  an  output  of  1572270  requisitions  the  Naval  Supply  Center 
would  have  an  output  of  1624129  requisitions.  Similarly,  for 
the  same  change  in  the  input  so  the  Naval  Supply  Center  the 
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output  of  PWC — Transportation  would  be  0.25  percent  greater 
than  the  baseline  output.  The  output  of  PWC — Transportation 
would  increase  from  1503492  gallons  of  fuel  to  1507196  gallons 
of  fuel.  Thus,  in  order  for  the  Naval  Supply  Center  to  output 
an  additional  51859  requisitions,  the  Naval  Supply  Center 
would  precipitate  the  use  of  an  additional  3704  gallons  of 
fuel  from  PWC — Transportation.  Further  the  Naval  Supply  Cen¬ 
ter's  change  in  output  would  cause  a  change  in  the  output  of 
Manpower — Civilian  of  0.07  percent  to  56929  manhours  of  an 
additional  38  manhours  of  civilian  labor. 

If  the  inputs  to  the  Naval  Supply  Center  were  ten  percent 
less  than  the  inputs  in  the  original  model,  the  percentage 
change  in  the  output  of  the  Naval  Supply  Center  would  be 
-3.29  percent.  This  would  imply  an  output  of  1520477  requi¬ 
sitions.  The  output  of  PWC — Transportation  would  be  1499794 
gallons  of  fuel  instead  of  1503492  gallons  of  fuel,  a  change 
of  -0.25  percent.  The  output  of  Manpower — Civilian  would  be 
56853  manhours  instead  of  56890  manhours,  a  change  of  -0.07 
percent.  The  above  figures  would  be  ranges  on  the  output  levels 
of  significant  sectors,  given  that  the  coefficients  for  the 
Naval  Supply  Center — San  Diego  were  within  plus  and  minus  ten 
percent.  These  ranges  would  be:  NSC--San  Diego  1520477 
to  1624129  requisitions;  PWC--Transportation  1499794  to 
1507196  gallons  of  fuel;  and  Manpower--Civilian  56853  to 
56929  manhours.  These  ranges  would  reflect  both  the  direct 
and  indirect  inputs  required  for  the  Naval  Supply  Center  to 
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operate  with  input-output  coefficients  within  plus  and  minus 
ten  percent  of  the  original  coefficients  of  the  model. 

The  administrator  of  the  Naval  Regional  Medical  Center 
might  wonder  what  services  would  be  needed  if  the  number  of 
shore-based  military  personnel  in  the  11th  Naval  District  was 
increased  or  decreased  by  ten  percent.  The  administrator  would 
only  have  to  go  to  Appendix  [B]  and  look  under  the  ten  percent 
column  of  Manpower — Active  Duty  and  read  off  the  percentage 
changes  for  the  two  sectors  of  the  Naval  Regional  Medical  Center 
and  see  that  there  would  be  0.19  percent  or  581  more  admissions 
for  a  total  of  310782  admissions  and  1.98  percent  or  34207 
more  outpatient  visits  for  a  total  of  1760873  outpatient 
visits.  The  minus  ten  percent  column  would  reveal  a  -0.19 
percent  change  or  an  output  of  309621  admissions  and  a  -1.98 
percent  change  or  an  output  of  1692562  outpatient  visits. 

Thus,  the  range  of  outputs  would  be:  Naval  Regional  Medical 
Center — Admissions  309621  to  310782  admissions;  and  Naval 
Regional  Medical  Center — Outpatients  1692562  to  1760873  out¬ 
patient  visits. 

Managers  of  the  various  sectors  of  this  model  might  use 
these  results  to  predict  future  requirements  for  their  sectors. 
It  was  pointed  out  above  that  predictions  could  be  made  as  to 
increased  requirements  based  on  increased  outputs  for  one  sec¬ 
tor.  Similarly,  decreased  requirements  and/or  resources  would 
result  from  a  decrease  in  the  input  to  a  sector  which  would 
necessitate  a  decrease  in  the  output  of  that  sector. 
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B .  RECOMMEND AT IONS 

Two  specific  areas  for  further  study  are  recommended. 

First,  the  model's  input-output  coefficient  matrix  needs 
to  be  reevaluated  to  see  if  all  the  zero  entries  are  still 
zero.  For  instance,  when  model  sectors  of  the  original  model 
formulation  were  split  to  form  other  sectors  were  the  inputs 
accurately  and  correctly  recorded?  One  would  think  that  per¬ 
haps  there  were  requisitions  submitted  to  the  Naval  Supply 
Center  in  support  of  outpatient  visits  to  the  Naval  Regional 
Medical  Center,  but  none  existed  in  the  model  as  formulated. 
Perhaps  it  was  not  possible  to  separate  the  supplies  that  went 
to  the  Naval  Regional  Medical  Center  because  of  the  manner  in 
which  accounting  was  done.  Nonetheless,  this  is  an  area  for 
further  study,  if  possible,  or  perhaps  a  means  of  tracking 
supplies  used  by  both  sectors  of  the  Naval  Regional  Medical 
Center  might  be  initiated  if  it  has  not  already  been  done. 

The  second  area  for  possible  further  study  is  all  sectors 
of  the  Public  Works  Center.  The  Public  Works  Center  was  the 
most  sensitive  sector  of  the  model  and  divisions  of  the  Public 
Works  Center  were  all  related  so  changes  in  the  coefficient 
matrix  here  could  be  crucial  to  the  model's  output.  Further 
work  in  this  area  could  prove  to  be  well  worth  the  effort  in 
terms  of  accurate  forecasts  by  the  model . 
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APPENDIX  A 

The  Technological  Coefficient  Matrix.  I-A 
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Tha  Final  Demand  Vector,  d 


Row  Sector 

1  NSC  -  SAN  CIEGC 

2  LBEACM  NAVAL  SHIPYARD 

3  INTER •  PAINT.  ACTIVITY 

4  NAVAIF  PEKRK  FACILITY 

5  PUBLIC  V.CFKS  CENTER 

6  PfcC  -  MAINTENANCE 

7  PMC  -  TRANSPQRTATI CN 

8  FV»C  -  UTILITIES  ENERGY 
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10  NAVRECNEDCEN  -  CUTFAT. 

11  NAVSTA  -  MILITARY  FERS. 
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17  NT C  -  SEPSCCLCCM 
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20  MANFCV.EP  -  CIVILIAN 


Demand 
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2687755. GO  MAN-FOURS 
6478818.00  MAN-HOURS 
194661.00  MAN-HOURS 
1074960.00  MAN-HOURS 
609420.00  GALLONS 
3751626.00  MILLION  0TUS 
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1499711.00  MAN-HOURS 
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759C4.00  PERSONNEL 
17856.00  PERSONNEL 
282C9.00  FERSCNNEL 


"Naval  SUPPLY  CEN?ER~ San  Diego 

Sector  -25%  -20  %  -15%  -10%  -5%  5%>  10%  15 %  20 %  25 % 

NSC  -  San  Diego  -8- 23  -6-58  -4-94  -J. 29  -1-65  1  65  3-30  4-95  6-60  8-25 
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MEAN  OF  X*S. 

VARIANCE  OF  X*S  (UNBIASED  ESTIMATE). 

STANOARO  DEVIATION  OF  X*S. 

T-STATISTIC  FOR  CONFIDENCE  INTERVAL  COMPUTATION. 
LONER  CONFIDENCE  LIMIT. 

UPPER  CONFICENCE  LIMIT. 

CONFIDENCE  INTERVAL  LENGTH. 

NLMBER  OF  OBSERVATIONS  MINUS  1. 

I  PLUS  1. 

£UMMY,  USED  IN  SORT  ROUTINE. 


START  PARAMETER  FOR  OUTPUT. 

STOP  PARAMETER  FCP  OUTPUT. 

N  CIV1DED  BY  2. 

N  CIVIDEO  BY  Z, 

N  CIVIDED  BY  2, 

N  DIVIDED  eY  4. 

N  CIVIDED  BY  4, 

N  CIVIDED  BY  4* 

LOWER  QUART IlE. 

MEOIAN. 

UPPER  GUARTILE. 

N  DIVIDED  BY  2, 

N  CIVIDEO  BY  2, 

RANGE  OF  X*$. 

INTERQUARTILE  RANGE. 

ftttiavSUEWtf!  cios5i“ciNT 

x*$. 

INTEGER  ROUNO  OF  K2. 

X*S. 

PERCENTILE  COEFFICIENT  OF  KURTOSIS. 


TIMES  2. 
PLUS  1. 

TIMES  2. 
PLUS  1. 


PLUS  N  DIVIDED  eY  4,  PLUS  1. 
PLUS  N  01 VI  CEO  EY  4,  PLUS  2. 


10  PERCENTILE  INDEX  OP 
90  PERCENTILE  INOEX  OF 


******************************* 


C 

C 

c 

c 

c 

c 

c 


CIMEKSION  A(  25*  25)*  AINV(25,  251 ,  AP(25,  25).  EASS(25),  0(25). 

1  CEL (25) *  OUTPUT (25.  50),  SECT(25,  24),  UNITSI25.  20), 

2  W  KAP  EA( 75  C ) ,  XRAY (25) 

PEAC  THE  INITIALIZING  PARAMETERS. 

READ  (4,  101)  IPRINT,  MSIZE,  NF,  NS 

I  WORK  «  (MSIZE  ♦*  2)  ♦  5  *  MSIZE 
MSIZE2  ■  2  ♦  MSIZE 

CALL  DRIVER  (A,  AINV,  Ap,  BASE.  0,  DEL,  OUTPUT,  SECT,  UNITS. 

1  VKAREA,  XRAY,  IPRINT,  IWORN,  MSIZE,  MSIZE2,  NF,  NS) 

WRITE  (NS,  102) 

STOP 

1?J  FORMAT  {*7«5  ///,  •  ALL  CONFUTATIONS  NOW  COMPLETE,  WITH  SUMMARY*, 

1  •  DATA  LISTED  ABOVE.  HAVE  A  NICE  DAY.*  T 

fieRCLTINE  CR1VER  (A,  AINV,  AP,  BASE,  0,  DEL,  OUTPUT,  SECT,  UNITS* 
1  WKAREA,  XRAY,  IPRINT,  IWORX,  MSIZE,  MSIZE2 ,  NF,  NS) 


97 


nnonn  r»o 


SECT (HSIZE*  241,  UNITSIMSIZE,  20),  WKAREA  IWORK) ,  XRAY(MSIZE) 


Zl  «  l.CE-6 


F  EAC  IK  THE  A  . 

CC  201  J  *  1,  M$ I ZE 

PE 40  (NF,  204)  (Ad, 
201  CCKTIKUE 


COEFFICIENT  MATRIX. 
MST" 


4  )*  I  -  1,  MSIZE) 


REAC 
RE  AO 


IK  THE  DEMAND  MATRIX. 

(NF,  2041  (C(I),  I  -1,  MSIZE) 


* 

* 


GIMEKSICNALITY  OF  5 
CF  24  IN  SECT (MSIZI 


4  COINCIDES  KITH  THE  DIMENSION 
,  24). 


* 

* 


REAC  EACH  SECTOR'S  IDENTIFICATION. 
CC  2£«_I_*_l,  MSIZE 


% 


REAC  (NF, 
2C2  CONTINUE 


2C5 )  (SECT( 1,  J),  J  ■  1,  24) 


REAC  EACH  SECTOR'S  UNITS. 

OC  2C3  1*1,  MSIZE 

REAC  (KF,  204 )  (UNITSd, 
2C3  CCKTIKUE 


J),  J  *  1,  MSIZE) 


E 


CALL  BMATRX  (A,  AF,  0,  SECT,  ZU,  MSIZE,  NF,  NS) 


WRITE  THE  HE AC  INC  ANO  PRINT  THE  A  AND  D  MATRICES. 

-  -  .  2(J7J 

208) 

MSIZE,  NF,  NS) 


n  r  4 

i  l 


( 

< 

CALL 


IPRINT 
I PRINT 
MPRINT 


.EC.  C 
.EC.  1 
(A,  0, 


NR 
)  HR 
SECT, 


HI  Hi: 


IPRINT, 


C 

C 


CALL  PARA  (A,  AINV,  AP,  EASE,  0,  CEL,  OUTPUT,  SECT,  UNITS 
1  WKAREA,  XRAY,  ZL,  IPRINT,  IWQRK,  MSIZE,  MSIZE2 ,  KF,  NS) 


PETUFK 


2C4  FCRMAT 
FCRMAT 


2C5 


2C6 
C  7 
08 


FORMAT 

FCRMAT 

FCRMAT 


(1F1C.C) 

MV 


('ITHf 
( ' 1THI 


CRICINAL  "A  AKO  0"  MATRICES:*) 
ORIGINAL  "A"  MATRIX:') 


!». 


Slercltike  para 

l  WKAREA,  XRAY, 


(A,  AINV,  AF,  BASE,  0,  DEL,  OUTPUT.  SECT,  UKITS 


i 


i«,  winy,  orjc,  u,  ucu,  uuirui,  ; _ 

ZL,  IPRINT,  I WORK,  MSIZE,  MSIZE2 ,  KF,  NS) 

CIMEN'IQN  AIMSIZE,  MSIZE),  AINVtMSIZE,  MSIZE).  AP(MSIZE,  MSIZE 

1  lAjElMSlZE),  C( MSIZE),  DEL  (MSI Z  E 2 1 ,  OUTPUTImS IZE ,  MSIZE2), 

2  ;ECT(M$ IZE,  24),  UNITS(MSIZE,  20),  WKAREA! IWORK ) ,  XRAY(MSIZE) 


CJMFjjTE.THE.aAS 


12  I 


00 


3C1 


2?} 

IF 

IF 


i’ 


M 

1, 


VECTOR  FCR  COMPARISON. 
IZE 


( 


CCKTINUE 
3C2  CONTINUE 
CALL  LINV2F 

“  1*5 

OC 


4  )  ■  A  ( 
I  .KE.  4 
I  .EC.  J 


J ) 


API  I, 
API  I, 


J  ) 
J) 


-AP(I,  J) 
1.0  -  AP(I, 


J ) 


(AP,  MSIZE,  MSIZE,  AINV,  1,  WKAREA ,  IEP) 


mm\  «  1,  MSIZE 
ASE(I)  -  C.O 


)  .LE.  ZL  ) 
A1NV( I ,  J) 


GO  TC 
*  D(J) 


303 


uuuouuu  u  u  u  u  u  u  u  u  u  uuuouuu 


CONTINUE 

CCKTINLE 

PEAC  THE  NUM6EF  Of 
PEAC  (NF*  3201  L II 
PEAC  (NF,  321  ) 


CELTA  POINTS  ANC  THE  DELTA  VALLES. 
(OEL(I),  I  -  1»  LIL) 


PPINT  TABLES  FCR  THE  THESIS  DRAFT. 

CALL  TABLE  (AP,  C.  SECT,  UNITS,  MSIZE,  NF,  NS) 

CC  RCN  AND  COLIMN  PARAMETRICS  FOR  EACH  OELTA  VALLE. 

cc  UM'f,:  IfatWtii' 

ICCLNT  «  0 

***************4****4^***4^****4i 
FARAMETPIC  ROW 

444*44*4*«*«***4**4>**«44*«**4*«* 

c;  310  I  -  1,  PSIZE 

CO  3C6  IL  ■  1,  HSIZE 

CO  3Cf  JL  ■  1,  MSI2E 

AF1IL,  JL)  «  A(1L,  JL) 

CONTINUE 

CONTINUE 

CO  207  J  *  1,  HSIZE 

IF  (  fiBS(AP(ltJI )  .LE.  ZL  )  GO  TO  307 
AP(  I  ,  J)  -  APU,  J)  ♦  DEL(LEE)  *  AP(I,  J) 

CONTINUE 

CO  209  11  «  1,  MSIZE 
CO  306  JJ  «  1,  MSIZE 

IF  (  II  .NE.  JJ  )  API 1 1  ,  JJ)  -  -AP (II,  JJ) 

IF  (  !l  .EQ.  JJ  )  AP(II,  JJ)  -  l.C  -  AP ( 1 1 «  JJ) 

CONTINUE 

CONTINUE 

CALL  LINV2F  (AP,  MSIZE,  MSIZE,  AINV,  1,  WKAPEA,  IER) 

CALL  FUL  (AINV,  0,  OUTPUT,  ZL,  ICOUNT,  MSIZE,  MSIZE2 ) 
CONTINUE 


FARAMTERIC  COLUMN 

CC  316  J  -  1,  FSIZE 
CO  212  11  -  1,  MSIZE 
CC  311  JL  -  1,  MSIZE 


cc  tthft  3l}’ •"*{!?,  JL, 

CONTINUE 


CONTINUE 
CO  3131 


GO  TO  312 
E)  *  AP  ( I,  J) 


99 


>0000000 


314 

315 


3U 


OC  315  II  -  It  MSIZE 

00  314  JJ  *  1,  MSIZE 
IF  I  II  .NE.  JJ  ) 
„  IF  III  .EQ.  JJ  ) 
CONTINUE 


API  1 1  f  JJ) 
APIII*  JJl 


CONTINUE 
CALL  LINV2F 


CALL  PUL  (AINVf 
CCNTINLE 


(AP,  HSIZEf  MSI ZEt  AINV« 

Of  OUTPUT,  ZL,  ICOUNT, 


•APIII.  JJ) 

1.0  -  APIII.  JJ) 


1,  WKAFEA,  I ER) 
MSIZE,  MSIZE2) 


IF  4gQ£LILEE) 

IF  4  OELILEEl 
♦  C.5) 

CALL  CCMPP  IAP, 

NS) 

IF  i  IPRINT  .NE.  0  .ANO.  IPRINT  .NE. 


•LT.  ZL  )  10ELTA  -  1FIXI (OEL ILEE)  *  100. C) 

•  CT.  ZL  )  IOELTA  «  IFIXI IOELILEE)  *  100. C) 

EASEf  OUTPUT,  SECT,  IOELTA,  MSIZE,  MSIZE2 ,  NF, 

1  )  60  TO  318 


317 


WRiIe  INS,  323)  MSIZE.  MSIZE2,  IOELTA 
CALL  LPRIN7  IOLTPUT,  SECT,  ICOLNT,  MSIZE,  NF,  NS,  l) 
{START  ■  1 

5uSpl?, 1 


IF  I  IPRINT 
CALL  XCGNF  IOU 
MSIZE2,  NF,  N'S  1 
ISTART 


GC  TO  316 
SECT,  XRAY,  IOELTA,  ISTART,  MSIZE, 


318 


C 

C 


316  CCNTINLE 
RETURN 


r  ,  r»  c 

MSIZE  ♦  1 
_  „  XCCNF  I OUTPUT , 
NSIIE2,  NF,  NS) 

^•0  TO 


jS  Pitt'.,,. 


SECT,  XRAY,  IOELTA, 
L1L  )  GC  TO  317 


MSIZE, 
ISTART,  MSIZE2,  MSIZE, 


320  FORMAT 


AT  II 
17  17 
AT  I  1 


15) 

7F10.0) 

PARAMETFICS  WITH 


321  FORM 

2*2  FORMAT  ( •  PARAMETFICS  WITH  A  DELTA  VALUE  OF*.  F6.2) 

323  FORMAT  ( •  ITHE  ',  12,  •  BY  12,  »  PARAMETRIZED  ROW'  * 

1  •  ANC  COLUMN  MATRIX  WITH  OSLTA  ■»,  14,  •  PERCENT.',  ///) 

SLBPCLTINE  XPRINT  IA,  0,  SECT,  IPRINT,  MSIZE,  NF,  NS) 

DIMENSION  AIMSIZE,  MSIZE),  CIMS1ZE),  SECTIMSIZE,  24) 

MMM4*4****4M*M*M4*4*****4** 

SLf ROUTINE  TO  WRITE  CUT  INTERMEDIATE  RESULTS  OF  "THESIS," 

WHICH  IS  A  PROGRAM  TO  CHECK  SENSITIVITY  IN  THE  I/C  MOOEL. 


« 

4 

IF  I 


IF  IPRINT  « 
IF  IPRINT  « 
IF  IPRINT  < 

IPRINT  .NE. 


0,  PRINT  BOTH 

1,  PRINT  CNLY 

2,  PRINT  CNLY 


THE  A  ANO  0  MATRICES. 
THE  A  MATRIX. 

THE  0  MATRIX. 


*  *  *  * 

* 

* 


C  .ANO.  IPRINT  .NE.  1  )  GO  TO  403 


\ 

C 

C 


SETUP  TO  PRINT  OLT  THE  MATRIX  PASSED  BY  COLUMNS  CF  10  ANC 
ffCWS  CF  MSIZE,  UNTIL  THE  MATRIX  tl  COMPLETELY  PRINTED  WITH 


4C1 


-  MSIZE,  UNTIL 
HEADERS. 

m'lUf  .GT.  MSIZE  ) 
WRITE  INS,  404 )  II,  I  - 

cc  in};  SH!1 


ISTOP  • 
ISTART, 


MSIZE 

ISTOP) 


it  It 

4C2  CCNTINL 


{All,  J),  J  •  ISTART,  ISTCPI 


100 


IF  <  ISTCF  .EC,  MSIZE) 
ISTAPT  «  ISTAPT  ♦  10 
ISTCF  «  ISTOP  4  1C 
WRITE  INS ,  40£) 

GC  TC  401 


€C  TO  403 


4C3  Tf  (  IPRINT  .NE.  C  .AND.  IPFINT  .NE. 

S  WRITE  CLT  THE  CEMAND  MATRIX  (MSI ZE  X 

WRITE  (NS,  405 I 


2  I 

1). 


RETURN 


C 

C 


C 

c 

c 

c 


4C4 


405 


4C6 
4CT 
4C8 
4CS 
4 1C 


CC  WRITE  7uk\  4li)E  ( SECT ( I ,  J),  J  ■ 
CONTINUE 


1,  24) ,  D<  I) 


RETURN 
FORMAT  (• 


*,  ////j  8X,  12,  5113) 


rurr »«  i  «  »  *  vwr  1 

FCRMAT  ('O',  1CF12.5) 
f?R?AT  <»1THE  CRIGINA 


.  ,  _  ORIGINAL  CEMANO  VECTOR:*,  //) 

FORMAT  4*0* ,2441,  P16.2) 

S LERCLT INE  BMATRX  (A,  B*  0,  SECT,  ZL,  MSIZE,  NP,  NS) 

CINENSION  AIMSIZE,  MSIZE),  E4MSIZE,  MSIZE),  OCMSIZE),  SECTIMSIZE* 
1  24) 

CC  5C2  I  *  1«  MSIZE 


5C1 


CC  501  J  -  1,  MSIZE 
fid,  J)  *  C.C 
IF  (  ABS ( A ( I »  J) 
efl,  j)  *  All,  J) 
CCNTIME 


.LT.  ZL 
*  0(J) 


)  GO  TO  501 


SC2  CCNTIME 


NF,  NS,  3) 


C 

C 

C 

C 

C 


5C3 


CALLELPPINT5(i!  SECT,  MSIZE,  MSIZE, 

PETLPN 

FORMAT  < • 1THE  INFLT  MATRIX,  e:»,  //) 

f LBR  CL  TI NE  UFFINT  (ARRAY,  SECT,  NCGLS,  MROWS,  NF,  NS,  ITYPE ) 

CIMENSICN  ARR AYIMPOWS,  NCOLS),  SECTfMROWS,  24) 

IF  (  ITYPE  .NE, 

ISTART  «  l 
ISTCF  *  5 
ICCLN  «  5 
GC  TC  601 


II  GC  TO  2 


ISTAPT 
ISTOP  > 
ICCLN  > 


•  1 
10 
10 


C 

c 


6C1  IF  <  ISTCP  .GT.  NCCLS  ) 


•Ji:  1 1 

60  5) 

CC  #Oi.I^S  MRCNS 


w :  i 

WPITE 


TYPE 

TYPE 

(NS, 


NR  IT 
WRIT 


ISTOP  « 


NCCLS 

6041 

610) 


:  i 


ISTART, 

ISTART, 


ISTCP) 

ISTOP) 


IF  ( 


ITYI 


.EG.  1 


1 


WRITE 


Avii,"j)i'j  *  t START,  IS Tup » 

IWE  ,lo.  2  I  WRITE  (NS,  608) 


Mil 


6C6) 


( SECT ( I ,  JI,  J 
(SECTU,  J),  J 


1, 

2, 


241, 

24), 


101 


nnnnftnon 


1  I  i,T«h|2iK  *<”' 

1  (APR AY (If  J),  J  ■  ISTART,  ISTOP) 

!  CONTINUE 

!cfc?uN  st  T0  “3 

UtCF  «  ISTOP  4  ICOUN 


( SECT (It  J),  J  «  It  24), 


tip  IT  f  (NS,  607) 
GC  TC  601 


6C2  FETIRA 

604  FORMAT  (•  •.  4CX*  12,  4120) 
At 5  FORMAT  j*0») 

606  FORMAT  (<0*,  24A1,  5F20.2) 

6C7  FORMAT  ('1*1 

608  FORMAT  (»0»,  24A1,  10F10.2) 

605  FORMAT  l*C*,  2441,  10F1C.0) 
610  FORMAT  (•  *,  3CX,  12,  9110) 


4CX*  12,  4120) 
24A1,  5F20.2) 


ItSpCLTINE  CON f  (>,  N,  NF,  NS) 

C  (PENSION  X(N ) 

********************************* 

SUBROUTINE  TO  CONFUTE  A  951  CONFIDENCE  INTERVAL,  GIVING  THE 
LOWER  ANO  UPPER  LIMITS,  FOR  THE  MEAN  CF  A  STRING  OF  NUMBERS, 


***** 


IF  (  N  ,LE.  1  )  RETURN 

JHalVo.o1 

21  •  1.0E-6 

cc  }IU  :  *6t!l  *  x c  i  * 

7C1  CONTINUE 

XEAR  «  TOTAL  /  FLOAT (N) 

VARX  •  0.0 

CC  (  X(I)  -  XEAR  )  ♦*  2 

7C2  CONTINUE 


VAPX  •  VARX  /  FLOAT (1C0F) 
SCEV  •  VARX  **  0.5 


kur  ,b i «  u 
COF  ,GT •  3 

o8|  rijr  I? 

CCF  JtI  40 


.ANO.  I  OOF  .LE. 
•ANC.  IDOF  .LE. 
.ANO.  IDOF  .LE. 
•  ANC.  I  OOF  .LE. 
.AND.  IOOF  .LE. 
.ANC.  IOOF  .LE. 
•ANO.  IOOF  .LE. 


»  I 

20  ! 
25  j 
30 

40  ) 


3 . 1  C 
2.36 
2.12 
2.09 

i:M 

2.C2 

2.CC 


XBAR  -  ( 
X8AR  ♦  I 
XU  "  XL 


T  *  SOEV  /  (FLOAT(N)  **  0.5 
T  *  SOEV  /  (FLOAT! N)  **  0.5 


ROUTINE  TC  SORT  TFE  ARRAY  ANO  PRINT  OUT  THE  LIST  CF  NUMBERS 
IN  THE  SORTED  ARRAY. 

NN1  •  N  -  1 

CC  104  I  -  1,  NMl 
IP1  •  I  4  1 


Lj 


on  on  n  oonnoon  n  ooooonon  onn 


C 

C 

c 


703 

7C4 


cc  xu, 

M W :  513 

X(JI  «  HCLC 
CCNT  IMJE 


1  GO  TC  703 


CCNTIM’E 
IF  (  N  .07 


100  )  GO  TO  706 


WFITE  (NS «  713) 


705 


J START  «  1 
IS7CP  »„10 


.GT.  N 
714) 


i 


HTOP  ■  N 

),  I  »  ISTART,  ISTOP  ) 
GO  TO  706 


7C6 


hT(PlSTCP 
WPI7E  (NS,  li*» 

IF  (  ISTOP  .EC.  N 
ISTART  *  I  START  ■» 

ISTCF  «  ISTOP  4  1C 
GC  TG  705 

hPITE  0LT  THE  STATISTICS  CCNPUTED  ABOVE. 

6FITE  (NS,  715)  N ,  IDOF,  TOTAL,  XEAR,  VARX,  SOEV,  T,  XL,  XU, 

4*4*4**************  *********** 

CONFUTE  CUARTILES  A  NO  COEFFICIENTS  OF  SKEWNESS  ANC  KURTOSIS. 
****************************** 


CL 


K  N  « 
NNP  1 

KI 


IS  N  EVEN  CR  CCC7 
NN  *  N  /  2 
NN  ♦  2 
«  NN  4  1 

COC  BRANCN  TC  708. 

FN  . LT.  N  )  GO  TO  7CE 

****** 

*  * 

♦  N  EVEN.  * 

*  * 
****** 

N  *  NN  /  2 
NN  •  IN  ♦  2 
N PI  •  IN  ♦  l 
S  N/2  EVEN  OR  Q0C7 
F  (  INN  .LT.  NN  ) 


BRANCH  TO  707  IF  ODD. 
GO  TO  707 


7C7 


N  EVEN  AKO  N/2  EVEN.  COMPUTE  QUARTILES. 
K1  «  NN  4  IN 

Ka  *  Kl  ♦  1 

Cl  -  (  X ( IN )  4  X(INPi)  )  /  2.0 

C2  >  (  X (NN )  4  X(NNP1)  )  /  2.0 

I  •  (  X(K1)  4  X ( K2 )  )  /  2.0 

C  TO  710 

COMPUTE  QUARTILES. 

)  /  2.0 


N  EVEN  ANC  N/2  CCC. 

K2  •  NN  4  INP1 

“  Vm\  ♦  X(NNP1) 

1  ic’W 


N  OOO.  * 
* 


103 

( 


Inpi  «  in  ♦  i 

IF  N il  000*  BP ANCF  TC  709. 

If  (  INN  .IT.  NN  )  GO  TO  709 


N  CCC  ANC  N/2  EVEN.  COMP 
K]  «  NN  4  IN  ♦  i 

KS  ■  K1  «  L 

Cl  »  (  X ( IN  1  4  X(INPl)  } 
C2  »  X(NNPl) 

C3  *  (  X  (K 1J  4  X(K2I  I  / 
GO  TC  710 


COMPUTE  QUARTILES. 


'  2.0 


N  CCC  AND  N/2  CCD. 
7C9  K2  »  NN  ♦  INP 1  4  1 
Cl  •  X(INPl) 

C2  »  X  I NNP 1 1 
C3  «  Xt  K2  ) 


COMPUTE  QUARTILES. 


CUARTILE  COMPUTATIONS  COMPLETE*  COMPUTE  SKEWNESS  AND  KURTOSTS 
COEFFICIENTS. 


7 1C  RANGE  »  X IN  1  -  X ( 1 1 
GRANGE  »  C2  -  Cl 
IF  (  AES (GRANGE)  .IT.  Zt 
CCCS  «  (C3  ♦  (-Z.C  *  Q2J 


CCCS  «  (C3  *  (-2.C  *  Q2)  ♦ 
IF  (  A ES (ORANGE )  .LT.  ZL  ) 
KFNC1  *  (0.1  ♦  FUCAT(N))  ♦ 
KPNC2  -  (0.9  *  FLCAT(N))  ♦ 
IF  (  AES(CPANGE)  .LT.  ZL  ) 


)  GO  TC  711 
4  Cl)  /  ORANGE 
)  QCOS  »  0.0 
4  C.5 
4  C.5 

)  GO  TC  712 


PCCK  «  (ORANGE  4  C.5)  /  (  XIKRND2)  -  X(KRNDl)  ) 
IF  I  AES (GRANGE)  .LT.  ZL  )  PCOK  *  0.0 

WFITE  CUT  THE  CONFUTED  STATISTICS  AND  STOP. 


WRITE  (NS*  716) 
I  FCCK 

RETURN 


X(N),  RANGE, 


GRANGE, 


713  FORMAT  <»  THE  SORTED  ARRAY  CF  NUMBERS  :  * ,  //  ) 

714  FORMAT  <*0*,  10F13.2  ) 

715  FORMAT  («  «,  ///,  •  STATISTICS:* »  //, 

1  *  N  «  *,  17,  Jt 

i :  55m  i’.  *», /, 

A  •  MEAN  «  *,  FIS. 2,  /, 

«  •  VAFX  -  *,  FIS. 2,  /, 

6  •  <CfV  ■  *t  FIS. 2,  /, 

7  •  ^-STAT  «  •*  F17.2,  /. 

€  »  LOWER  CCNF ICENCE  LIMIT  «  *,  F19.2,  /, 

S  «  UPPER  CCNF  ICENCE  LIMIT  ■  »f  F19.2,  //, 

«  •  CONFIDENCE  INTERVAL  LENGTH  *  *,  F15.2  ) 


716  FCFMAT 
1  •  MAX 


vrrcn  vunriucnuc  Linn  -  -*  r/f 

CONFIDENCE  INTERVAL  LENGTH  *  *,  F15.2  ) 

MAT  (»0MINIMUM  *  *»  F20.2*  /, 

MAXIMUM  -  ».  F2C.2#  /, 

RANGE  •  •«  922*2 1  /, 

LCWfR  CUARTILE  »  *,  F13.2,  /, 

MEDIAN  -  *,  F21.2,  /, 

YKfK-Stftit!  /. 

CUARTILE  COEFFICIENT  OF  SKEWNESS  ■  •*  F9.4,  / 
PERCENTILE  COEFFICIENT  OF  KURTCSIS  «  «,  F7.4, 


SUBROUTINE  XCCNF  (ARRAY*  SECT*  XRAY,  IDELTA*  I  START,  ISTOP, 
l  MSI Z E *  MSIZE2*  NF*  NS) 

CIMENSICK  ARRAY (MS  IZE*  MSIZE2),  SECTIMSIZE,  24),  XRAY(MSIZE) 
IF  C  ISTAPT  .LT.  NSIZE  )  WRITE  (NS,  803)  IDELTA 


IF  (  J START  .GT.  ►SIZE  )  WRITE  (NS,  804)  IDELTA 


CC  802  I  > 
JJ  «  0 


MST2S 


CO  SCI  J  •  I  START ,  ISTGP 
JJ  «  JJ  ♦  1 

XRAY(JJ)  A  ARRAY! 1,  J) 

8C1  CONTINUE 

IF  (  I  .GT.  l  I  WRITE  (NS,  80S) 

WRITE  (NS,  806 )  (SECT! I,  JJ),  JJ  -  1,  24) 

CALL  CCKF  (XRAY,  MSIZE,’ NF,  NS) 

8C2  CONTINUE 

RETURN 

8C3  FORMAT  ( *1$TAT 1ST  ICS  ON  PARAMETRI2E0  ROWS  WITH  CELTA  «',  14, 

1  •  PERCENT.*,  /✓/ 1 

8C4  FORMAT  ( • 1 STAT I  ST  ICS  ON  PARAMETRIZED  COLUMNS  WITH  CELTA  »*,  14, 
1  «  PEPCENT.*,  ///) 

£05  FORMAT  4*1*1 


8C6  FORMAT  (*  * 


DEMAND  FOR - •  ,  24A1 ,  *  * 


///) 


SLERCLTINE  CCMFR  (AP,  BASE,  OUTPUT,  SECT,  IDELTA,  MSIZE, 

1  MSIZE2,  NF,  NS) 

CIMENSICN  AP ( MS IZ E,  MSIZE),  BASE(MSIZE),  OUTPUT! MSIZE,  MSIZE2), 

1  SECT (MS  IZE,  24),  AVGROW42S),  AVGCCL! 25) ,  PCTR(25),  PCTC(25J 

WRITE  (NS,  910)  ICELTA 

CC  9C2  I  «  1,  MSIZE 
TOTAL  *  C.C 

CC  901  J  »  1,  MSIZE 

TOTAL  -  TOTAL  +  OUTPUTd,  J) 

90 1  CONTINUE 

AVGROWd)  *  TCTAL  /  FLOAT (MSI ZE) 

TCTAL  *  0.0 
ISTART  ■  MSIZE  •*  1 

CO  902  J  *  ISTART,  MSIZE2 

„  TCTAL  »  TOTAL  ♦  OUTPUTd,  J) 

902  CONTINUE 

AVOCCL(I)  »  TCTAL  /  FLOAT(MSIZE) 

9C3  CCNTINLE 

WRITE  (NS,  911) 

CC  9C4  I  >  1,  MSIZE 

PCTRd)  »  (  (AVC-ROWd)  /  BASE(I))  *  100.0  )  -  1C0.0 
FCTCdJ  -  (  (AVGCOL(I)  /  EASE!  I ) )  4  100. 0  )  -  100.0 
WRITE  (NS,  912)  (SECT!  I,  J),  J  ■  1,  24),  6ASE(I),  AVGRCW(I), 
1  PC  TP  ( I ) ,  AVGCCL(I),  PCTC(I) 

904  CCNTINLE 

WRITE  (NS,  913)  ICELTA 

CC  90S  I  ■  1,  MSIZE 
I]  •  I  ♦  MSIZE 

PCTR(j)  •  (  (OUTPUT! 1,  I)  -  BASE(I)  )  !  BASE(I)  )  *  100. C 
PCTC(I)  •  (  (OUTPUTd,  II)  -  8ASE(I)  )  /  8ASE(I)  )  *  100.0 
WRITE  (NS,  914)  (SECT  (  I ,  J),  J  »  1,  24),  EASEd), 

1  OUTPUT  ( I,  I),  FCTRd),  OUTPUTd,  II),  PCTC(I) 

5C5  CCNTINLE 

CC  SC9  J  •  1,  MSIZE 
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JJ  •  J  ♦  MSIZE 

WRITE  (NS,  913)  IOELTA 

WRITE  (NS,  915 )  J,  J ,  J, 


JJ 


CC 


*??R( 1 


C 

C 


506 


PCTR(I)  •  (  f OUTPUT  <  I ,  J )  -  EASE(I)  )  /  BASE(I) 
IF  (  IDELTA  .EQ.  10  )  AF(I,  J)  »  PCTR(I) 
PCTC(I)  «  (  (OUTPUT (I j  JJ)  -  BASEd 
WRITE  (NS,  514)  ( 

OUTPUT!  I«  J),  PCTRi 
CONTINUE 


)  * 


T (I  ,  JJ)  -  BASE(I)  )  /  EASE(I)  )  * 
(SECTd,  K).  K  -  1,  24)  ,  BASE(I), 
(I) ,  OUTPUT ( I ,  JJ),  PCTC(I) 


1C0.0 

100.0 


TCTP 

TCTC 


0.0 

c.c 


c 

c 
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CC  907  K  «  1,  MSIZE 

TCTP  »  TCTP  ♦  PCTR(K) 

TOTC  «  TCTC  ♦  PCTC(K) 

PMEANR  «  TCTR  /  FLOAT (K) 

PMEANC  •  TCTC  /  FLOAT (K) 

CONTINUE 

WRITE  (NS,  5161  PMEANR 

CC  506  K  *  1.  MSIZE 

IF  <  ABS(PCTR(KI)  .GT.  A8S I PMEANR)  ) 


5CE 


(SECT (K ,  L),  L 
CONTINUE 


-  1,  24),  PCTR(K),  PMEANR 


WRITE  (NS,  917)  K, 


C 

C 


909  CONTINUE 

IF  (  ICELTA 
IF  (  ICELTA 
1  NF,  NS,  2) 


.EC.  1C  I  WRITE  (NS,  918) 
.EC.  10  )  CALL  UPRINT  (AP, 


SECT,  MSIZE,  MSI ZE , 


510 

511 

912 

913 

914 

515 

516 

517 

5ie 


PETLPN 

FCPMAT  (»1PCW  ANC  COLUMN  MEANS  FOP  EACH  SECTOR  WITH  DELTA  »', 

1  14,  •  PERCENT.*,  ///) 

FORMAT  (*  *,  'SECTOR*.  27X,  ‘BASE*,  19X,  'ROW  MEAN*,  6X, 

1  *t  OIF  (ROW)*,  11X,  *CQL  MEAN*.  6X,  'X  OIF  (COL)*) 
lFOPFAT5,.C.,  24AI,  5X,  F12.2,  12X,  F12.2,  7X,  F*.2,  13X,  F12.2, 

1RCFPAT5d  1EFFECT  CF  ROW  PARAMETRICS  ON  EACH  SECTOR  WITH*, 

1  •  DELTA  14,  '  PERCENT.',  ///,  '  SECTOR*,  27X, 

2  'EASE',  18X,  'ROW  OUTPUT',  5X,  «X  DIF  (SEC)1, 

3  n>,  «CCU  OUTPUT',  5X,  *1  CIF  (CCL)') 

FORMAT  ('O',  24A1,  5X,  F12.2,  12X,  F12.2,  7X,  F7.2,  11X,  c12.2, 

1  7X,  F7.2) 

FORMAT  (  'GROW * ,  12,  '  -  COLUMN',  13,  *  OF  OUTPUT', 

1  •  ♦  *  ♦  *  COL',  12,  •  -  COLUMN',  13,  •  OF  OUTPUT.)') 

FCPMT  CC,  //,  «  SUMMARY  FOR  PARAMETRICS  CN  PCWS:  MEAN  »  ', 

1  FE.4 ,  »  X.') 

FORMAT  CCPERCENTAGE  CHANGE  FOR  SECTOR  ',  12,  ',',  IX,  24A1,  /, 

1  •  CF  ',  F6.2 ,  '*«, 

CPEATEP  IN  AeSCLUTE  MAGNITUDE  THAN  THE  MEAN  PERCENTAGE  CFANGE* , 


FCRMaI.  CCPERCENTAGE  CHANGE  FOR  SECTOR  ',  12,  ',', 

« 

2  •  OF  *,  F6.2,  '!.*) 

FORMA  1  C1MATPIX  CF  PERCENTAGE  CHANGES  WITH  DELTA  EQUAL  TC  1C*:', 
1  //) 


C 

c 

c 


ENC 

SUBROUTINE  MLl  (AINV,  0,  OUTPUT,  ZL,  ICOUNT,  MSIZE,  MSIZE2) 

CIMENSION  AINV (MSIZE,  MSIZE),  D(MSIZE),  OUTPUT (MS  IZE ♦  MSIZE2) 

ICCUNT  «  ICOUNT  ♦  1 

CC  1002  I  »  1,  MSIZE 

OUTPUT (I ,  ICCUNT)  «  0.0 

CO  ICC  1  J  *  1,  MSIZE 

IF  (  ABS( A1NV( I ,  J))  «UE.  ZL  I 


GO  TO  10C1 
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1 


OUTPITII,  ICCUNTI 
.1001  CONTINUE 

10C2  CCNTIALE 


»  OLTPUT  ( I  *  I  COUNT)  ♦  AINVd,  J)  « 


C 

C 


PETIPA 
E  AD 

SlESOUTINE  TABLE  (AP,  0,  SECT,  UNITS,  MSIZE,  NF,  NS) 
iCIjjEAJJCAjAPCPJIZE,  MSIZE),  D(MSIZE),  SECTIMSIZE,  24), 

♦•FITE  TABLES  CF  "(I-A)"  ANC  "D"  MATRICES. 

♦.FITE  (NS,  11C3) 


I  .EC.  21  ) 


C 

c 

c 

c 


cc  i?',1  i  :d:  iv% _ 

HOI  CCHlSJI  <NSl  11041  '*  ,4F<1'  J1’ 
♦.BITE  (NS,  UC2) 


WRITE  (NS,  11C3) 
1,  MSIZEI 


CC  1 1 C 2  I  *  1,  MSIZE 

WRITE  (  NS,  11C2)  ( S£CT(  I,  J),  J  ■  1,  24),  0(1), 
1  (UAITSd,  J),  J  -  1,  20) 

11C2  CCA.TIACE 

FETLPA 


1103  FORMAT  («1»,  /////) 

11C4  FCFMAT  («  •,  26X,  »ROW(»,  12,  *)*,  7(/,  4F14.5 )) 
11C5  FCPMA7  (•€',  4X,  Z4A1,  2X,  fll.2,  IX,  20A1  ) 

EAC 


C(  J) 
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